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Abstract 



After an introduction to the problem of cosmological structure formation, we develop gauge invari- 
ant cosmological perturbation theory. We derive the first order perturbation equations of Einstein's 
equations and energy momentum "conservation" . Furthermore, the perturbations of Liouville's equa- 
tion for collisionless particles and Boltzmann's equation for Compton scattering arc worked out. We 
fully discuss the propagation of photons in a perturbed Friedmann universe, calculating the Sachs- 
Wolfe effect and light deflection. The perturbation equations are extended to accommodate also 
perturbations induced by seeds. 

With these general results we discuss some of the main aspects of the texture model for the 
formation of large scale structure in the Universe (galaxies, clusters, sheets, voids). In this model, 
perturbations in the dark matter are induced by texture seeds. The gravitational effects of a spher- 
ically symmetric collapsing texture on dark matter, baryonic matter and photons are calculated in 
first order perturbation theory. We study the characteristic signature of the microwave background 
fluctuations induced in this scenario and compare it with the COBE observations. 
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Introduction 



Within standard cosmology, the formation of large scale structure remains one of the biggest un- 
solved problems, despite of great efforts. The most natural scenario, where structure forms by 
growth of adiabatic perturbations in a baryon dominated universe is clearly ruled out. Other exten- 
sively worked out scenarios like isocurvature baryons or baryons and collisionless matter (cold dark 



matter or hot dark matter) face severe difficulties [ Seller and Huchra, 1989 , Maddox et al., 1990 



paunders et al., 1991| , |Gouda et al., 199"l| (for a short review see [ Peebles and Silk, 1990| ). 

On the other hand, there is the interesting possibility to induce perturbations in the baryonic and 
dark matter by seeds. Seeds are an inhomogeneously distributed form of energy which contributes 
only a small fraction of the total energy density of the universe. Thus, linear perturbation theory 
can be used to calculate the induced fluctuations and their time evolution. Gauge-invariant linear 
perturbation theory [ Bardeen, 1980| [ is superior to gauge dependent methods, since it is not plagued 
by gauge modes, and it leads in all known cases to the simplest systems of equations. 

Examples of seeds are primordial black holes, boson stars, a first generation of stars, cosmic 
strings, global monopoles, or global texture. These last three are especially appealing since they can 
originate in a natural way from phase transitions in the early universe [[Kibble, 198q] . If they indeed 
play an important role in structure formation they would link the smallest microscopical scales never 
probed by particle accelerators (down to W~^^GeV~^) and the largest structures (up to lOOMpc 
and maybe more)! Another attractive feature of topological defects is that the gravitational effects 
of each class of them (cosmic strings are tti defects, monopoles are tt2 defects and texture are tts 
defects) are quite insensitive to the detailed symmetry breaking mechanism and only depend on the 
symmetry breaking scale t], which we cast in the dimensionless quantity e = IGvrGry^. This quantity 
determines the amplitude but not the shape or the time evolution of the perturbations. 

Here, we mainly consider the defect global texture which was first proposed by Turok [1989] 
as a seed for large scale structure. Several subsequent investigations of this scenario gave promising 
results: The spatial and angular correlation functions, the large scale velocity fields and other sta- 
tistical quantities obtained by numerical simulations agree roughly with observations (see Pen et al. 
[1993], Gooding et al. [1992], Spergel et al. [1991] and references therein). 

My objective in this text is to fully develop gauge invariant linear perturbation theory to treat 
models with seeds. I then want to show in some detail, how all the linear perturbation theory 
aspects of a scenario of large scale structure can be investigated with these tools. As an example, 
we discuss the texture scenario. I choose this scenario not because I think it is the solution of the 
problem of cosmological structure formation. But it is the simplest worked out scenario, where 
initial fluctuations are induced by topological defects of a symmetry breaking phase transition, and 
I believe that this class of models deserves thorough investigation as an alternative to models with 
initial perturbations from inflation. 

In the first chapter we set a frame for this review with a non-technical overview of the problem 
of cosmological structure formation. We also very briefly discuss some of the presently considered 
scenarios. A reader already familiar with the problem of structure formation may skip it. In the 
second chapter gauge-invariant cosmological perturbation theory [ Hawking, 1966| , Bardeen, 1980 



[Kodama and Sasaki, 1984 , purrer and Straumann, 1988| , [Mukhanov et al., 1992| ] is presented in a 
(hopefully) pedagogical fashion. Although, some familiarity with general relativity and basic concepts 
of differential geometry are needed to follow the derivations in this chapter. We extensively discuss 
scalar, vector and tensor perturbations for fluids and for collisionless matter. Furthermore, general 
formulae for the deflection of light in a perturbed Friedmann universe are derived. This part is new. 
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Finally, seeds as initial perturbations are introduced [Durrer, 1990]. Here it is assumed that non- 
gravitational interactions of the seeds with the surrounding matter are negligible. For topological 
defects, this is certainly a good approximation soon after the phase transition. 

In Chapter 3, some simple but important applications of cosmological perturbation theory are 
described. We first discuss the ideal fluid [ Bardeen, 198C , Kodama and Sasaki, 1984 1. Then, a gauge 
invariant form of Boltzmann's equation for Compton scattering is derived. We study it in the limit of 
many collisions to obtain an approximation for the damping of cosmic microwave background (CMB) 
fluctuations by photon diffusion in a reionized universe. We then list all mechanisms proposed to 
induce anisotropics in the cosmic microwave background. Finally, two applications of light deflection 
are discussed: Lensing by global monopoles [ Barriola and Vilenkin, 1989 1 and light deflection due to a 
passing gravitational wave. This effect is new and might be an alternative way to detect gravitational 
waves of very far away sources. 

In Chapter 4 we review the concept of texture defects and present the exact flat space solution 
found by Turok and Spergel [1990]. We then calculate the gravitational potentials induced by this 
solution and apply the formalism developed in Chapter 2 to derive photon redshift, light deflection 
and perturbations of baryons and dark matter. All results can be obtained analytically and provide 
a nice application of gauge invariant perturbation theory. They are, however, cosmologically relevant 
only for scales substantially beyond horizon scale. Most of the results of this chapter were originally 
derived by Durrer [1990] and Durrer et al. [1992a] , but the calculation of light deflection is new and 
the work on collisionless particles is presented in a much simpler way and a physically more sensible 
limit is performed. 

We conclude with a chapter on textures in expanding space. Here we present a method to 
calculate the induced cosmic microwave background fluctuations by statistically distributing indi- 
vidual textures which are modeled spherically symmetric. The detailed numerical results of these 
investigations will be published elsewhere ]purrer et al., 1993 ], 

In Appendix A the basics of 3+1 formalism of general relativity are outlined. In Appendix B, we 
provide a glossary of the variables used in this review. 
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Chapter 1 



The Problem of Large Scale Structure 
Formation 



1.1 The Standard Cosmological Model 

We begin with some facts on the standard model of cosmology. Extensive treatments of this subject 
can be found in the books by Tolman [1934], Bondi [1960], Sciama [1971], Peebles [1971], Weinberg 
[1971], Zel'dovich and Novikov [1983], Borner [1988], Kolb and Turner [1990], Peebles [1993]. 

The assumption that the universe is homogeneous and isotropic on very large scales, leads to a 
very special class of solutions to the gravitational field equations, the Friedmann-Lemaitre universes. 
Accepting furthermore the cosmological origin of quasars and/or the cosmic background radiation, 
one can reject the possibility of a so called bouncing solution and show that our universe necessarily 



started in a big bang [ Ehlers and Rindler, 198£ ]. Due to homogeneity and isotropy, the energy 



momentum tensor is described by the total energy density p and pressure p which are functions of 
time only. The metric is of the form 

ds^ = a^i-df +-fijdx'dx^) , 

where 7 is a 3-d metric of constant curvature k. The simply connected 3-spaces corresponding to 
7 depend on the sign of the curvature: For /c > 0, the three-space is a sphere, for A; < a pseudo 
sphere and for /c = flat, Euclidean space. (Note, however, that the metric 7 cannot decide on the 
topological structure of 3-space. For A; = 0, e.g., it may well be, that three space is topologically 
equivalent to R^/Z^, i.e., a torus with finite volume. The often stated phrase that for /c < 3-space 
is infinite is thus wrong.) 

The time dependent function a is the scale factor and the physical time (proper time) r of an 
observer at rest is given by dr = adt. The time coordinate t is called conformal time. 

Einstein's equations imply the Friedmann equation which determines the scale factor as a function 
of the density: 

. 2 , A„2 



- +k = —Gpa^ + Aay3 . (1.1) 
a J 3 

Here G is Newton's gravitational constant, and A is the famous cosmological constant which has 



been resurrected several times in the past [ Weinberg, 198£ ] 
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Important quantities in Friedmann cosmology are the Hubble parameter, H and the density 
parameter Q, which are defined by 

H = ale? and = p/p^ , with p^ = 3H^/{8itG) 

. Here pc is the critical density, the density of a universe with A; = A = 0. The present value of the 
Hubble parameter, Hq, is usually parameterized in the form Hq = h x lOOkm/s/Mpc. Observations 
limit the value of Hq in the range 

0.4 < /i < 1 and 
0.05 < Oo < 2 . 

From estimates of the mass to light ratio one obtains a value for the amount of luminous matter, 
^lum ~ 0.007 and from velocity measurements one estimates the amount of 'clustered matter' to be 
^dyn ~ 0.1. These values result from many difficult observations and are correspondingly uncertain 
[Kolb and Turner, 199C]. Unfortunately, the measurements which lead to them can not provide 



information about a dark, poorly clustered matter contribution and therefore do not yield an upper 
limit for the density parameter. The upper bound comes from a comparison of lower limits to the 
age of the universe with the Hubble parameter. 
Using the energy "conservation" equation 

p = -3-{p + p), (1.2) 
a 

one finds that for p> 0, p grows at least like a~'^ for a ^ (i.e. approaching the big bang). Therefore, 
at an early enough epoch the density term in ( |1.1D always dominates over the curvature term and 
the cosmological constant, and becomes arbitrarily close to 1. In other words, the evolution of a 
Friedmann universe always "starts" very close to the unstable fixpoint Q = 1. 

On the other hand, ( |1.1| ) shows, if the cosmological constant A once dominates the expansion of 
the universe, the curvature term k and the energy density p become less and less important. The 
universe then expands approximately exponentially in terms of physical time r 



a(T) = a(rj) exp(Y^A/3 / dr') 

J To 



This rapid expansion can only be stopped if some or all of the vacuum energy inherent in the cos- 
mological constant is radiated into particles. Such an intermediate phase of A dominated expansion 
is the basic idea of most scenarios of inflation. 

Today, the cosmological constant is severely limited from observations: 

This is very much smaller than the values expected from particle physics, but might still be enough 
to influence the expansion of the universe and structure formation substantially pioltzman, 1989 1. 



It can, e.g., lead to a 'loitering period' during which the universe is nearly non-expanding 

[Durrer and Kovner, 1990 1, and thus fluctuations grow exponentially. Since such a small value for A 



requires an enormous fine tuning, we set A = in the following. 

For A = 0, one has > 1 in a closed universe (3-sphere) and $7 < 1 in a negatively curved (open) 
universe. 
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Later, we also use the following consequences of (IJ.) and ( |1.2| ) (setting w = p/ p and = p/ p)'. 

-\[{3w-i)C-f + {3w + m (1.3) 

w = 3{d,/a){w-ci){l + w) . (1.4) 



a/a 



Like the standard model of particle physics, the standard model of cosmology has many impressive 
successes. The most important ones are: 

• Uniform Hubble expansion (Detected by Wirtz [1918], Shpher [1920] and Hubble [1927]). 

• The prediction of Gamov [1946] and detection by Penzias and Wilson [1965] of the cosmic 
microwave background radiation, its 'perfect' black body spectrum [ Mather et al., 1990|] and 
its extraordinary uniformity (see Figs. 1 and 2). 

• The abundance of the light elements {H, ^H, ^H, ^He, '^He, Li) can be calculated, 
and the comparison with observational estimates predicts 0.02 <VtB< 0.1 which is consistent 
with direct determinations of VLq. This calculations were originally performed by Alpher et al. 
[1948], Wagoner et al. [1967] and by many others later. A comprehensive review is Boesgaard 
and Steigman [1985]. For more recent developments, see e.g. Kurki-Suonio et al. [1988] and 
Walker et al. [1991]. 

However, many questions are left open. Most of them can be cast in terms of very improbable initial 
conditions: 

• The question of the cosmological constant. A: Why is it so small? i.e., so much smaller than 



typical vacuum energies arising from particle physics [ Weinberg, 1989 ] 



• The flatness / oldness / entropy problem: why is the universe so old, t » tpu and still 0(57) = 1? 
Here, tpi is the Planck time, tpi = {fiG/c'Y/'^ = 5.39 x lO'^^^sec . 

• The horizon problem: how can different patches of the universe have been at the same tem- 
perature long before they ever where in causal contact? 

So far, the most successful approach to answer the second and third questions is based on the 
hypothesis of an inflationary phase during which expansion is dominated by a 'fluid' with neg- 
ative effective pressure, p < — p/3 (e.g. a cosmological term, for which p = —p) and physical 
distances thus grow faster than the size of the effective particle horizon, (d/a^)~^. I do not ex- 
plain how inflation answers these questions, but refer the reader to some important publications 



on the subject: jGuth, 198l| , [Albrecht and Steinhardt, 1982| , |Linde, 198^ |Linde, 1984 |Linde, 1990 



La and Steinhardt, 19891 Plive, 1990 |. 

These open problems provide some of the very few and important hints to new fundamental 
physics beyond the standard models of cosmology and of particle physics. 



1.2 Structure Formation 

Besides the fundamental issues, related to the very early universe and unknown basic physics, there 
is the important problem of structure formation which, to some extent, should be solvable within 
the standard model. This is the problem I want to address here: 
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How did cosmological structures like galaxies, quasars, clusters, voids, sheets... form? 
From difficult and expensive determinations of the three-dimensional distributions of optical and 
infrared galaxies on scales up to 150/i~^Mpc [De Lapparent et al., 1986, De Lapparent et al., 1988, 



Peller and Huchra, 1989, Broadhurst et al., 1990 , paunders et al., 1991 , Loveday et al., 1992a |, 



[Fisher et al., 1993 1, we know that galaxies are arranged in sheets which surround seemingly empty 



voids of sizes up to 50h ^Mpc. Observations also show that the distribution of clusters of galaxies 



Bahcall and Soneira, 1983, Bahcall, 1989, West and van den Bergh, 1991] is inhomogeneous with 



6nc 

> 1 on scales up to 25 Mpc. 

In addition, it is important to observe that not all galaxies are young. There are quasars with 
redshift up to z ~ 5 and galaxies with z ~ 3. On the other hand, there are indications that the 
galaxy luminosity function is still evolving considerably in the recent past, i.e., that galaxy formation 
still continues |]Loveday et al., 1992b |. 




Prom a naive Newtonian point of view one might say: "Gravity is an unstable interaction. Once 
the slightest perturbations (e.g. thermal) are present they start growing and eventually form all these 
observed structures. The details of this process might be complicated but there seems to be no basic 
difficulty. " 

This is true in a static space where small density fluctuations on large enough scales grow ex- 
ponentially. But, as we shall see in Chapter 2, in an expanding space, this growth is reduced to a 
power law: 

for pressureless matter, dust 
for radiation, relativistic particles. 

In a radiation dominated universe, pressure prohibits any substantial growth of density perturba- 
tions. Once the universe is matter dominated 5p/p grows proportionally to the scale factor, but the 
gravitational potential 

^' oc y {6p/r)r^dr = J Sprdr 

remains constant. This led Lifshitz [1946], who first investigated the general relativistic theory of 
cosmological density fluctuations, to the statement: "We can apparently conclude that gravitational 
instability is not the source of condensation of matter into separate nebulae". Only 20 years later 
Novikov [1965] pointed out that Lifshitz was not quite right. But to some extent Lifshitz's point 
remains valid: Gravity cannot be the whole story. Imagine an Q = 1 universe which was matter 
dominated after z ~ 2 x 10^. Since density perturbations in cold matter evolve like l/{z + 1), they 
have grown by less than 10^. But today there are galaxies with masses of M « lO^^M© ~ lO^^mN ■ 
Statistical inhomogeneities of these N nucleons have an amplitude of 1/VN « 10~^^, which would 
have grown to less than 10~^^ <^ 1 today. Therefore, there must have been some non-thermal initial 
fluctuations present with an amplitude on the order of 10^^ — 10~^. We also know that the amplitude 
of these initial perturbations did not exceed this amount because of the high degree of isotropy of 
the cosmic microwave background radiation (CMB). In fact, only recently CMB anisotropies have 
been found in the COBE experiment [ Wright et al., 1992| , Smoot et al., 199^ ]. They are on the level 



AT 

— 10"^ or (5/3(^ad)Mrad) ~ 4 X 10"^ 
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on all angular scales larger than 10° and compatible with a scale invariant Harrison-Zel'dovich 
spectrum. On smaller angular scales other experiments led to limits, AT/T < a fewxlO^^ (see 
Fig. 2). We disregard the dipole anisotropy which is due to our motion with about 600km/s with 
respect to the microwave background radiation. 

There are scenarios of cosmological structure formation which are not based on the gravitational 
instability picture. The most recent of them is based on explosions (of supernovae or super conducting 
cosmic strings) fOstriker and Cowie, 1981 ] which are supposed to sweep away baryons, producing 
shock fronts in the form of sheets surrounding roughly spherical voids. Usually such scenarios face 
enormous difficulties in producing the energy required to account for the large scale observations and 
to satisfy at the same time the limits for perturbations of the CMB on small angular scales. None 
of them has thus been worked out in detail. In fact, the explosion scenario also needs gravitational 
instability to amplify perturbations by a significant factor |]McKee and Ostriker, 1988 ] . In this paper, 
we concentrate on scenarios which rely on gravitational instability. The fact, that the anisotropies 
measured by COBE just coincide with the amount of growth necessary to form structures today is 
taken as a hint that the gravitational instability picture may be correct. 



1.3 The General Strategy 

We now want to outline in some generality the ingredients that go into a model of structure formation 
which is based on the gravitational instability picture. 



1.3.1 Initial fluctuations 

We saw in the last section that small density fluctuations in a Friedmann universe may have grown 
(by gravitational instability) by about a factor of 2 x 10^ during the era of matter domination. There- 
fore, a complete scenario of structure formation must lead to initial matter density fluctuations with 
amplitudes on the order of 10~^. Two possibilities to obtain these initial fluctuations are primarily 
investigated. 

A) Initial perturbations produced during inflation: Here it is assumed that density 
fluctuations are generated during an inflationary phase |Guth, 1981, Albrecht and Steinhardt, 1982 



Linde, 1982 , Linde, 1983 , Linde, 1984| , La and Steinhardt, 1989 , Linde, 1990f| from initial quantum 



fluctuations of scalar flelds [Bunch and Davics, 1978, Hawking, 1982, Starobinsky, 1982], 
[ Guth and Pi, 198^ , Fischler et al., 1985 ]. Due to the nature of quantum fluctuations, the distribu- 
tion of the amplitudes of these initial perturbations is usually Gaussian. 

There are several difi^erent more or less convincing models of inflation. One divides them into: 
Standard (or old) inflation ] | Guth, 1981 ], new inflation ]Albrecht and Steinhardt, 1982, Linde, 1982 



chaotic inflation ]lLinde, 1983 ] and the most recently proposed possibilities of extended and hyper- 
extended inflation ]|La and Steinhardt, 1989[] . Reviews on the subject are found in Linde [1984], 
Linde [1990], Olive [1990] and Steinhardt [1993]. All these models of inflation differ substantially, 
but the mechanism to produce initial fluctuations is basically the same: 

The scales / of interest {I < Inito) ~ 3000/i^^Mpc) are smaller than the effective particle horizon at 
the beginning of the inflationary phase. The scalar field that drives inflation therefore experiences 
quantum fluctuations on all these scales. Their amplitude for a minimally coupled scalar fleld 



m 
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de Sitter universeFI within a volume V can be calculated [Bunch and Davies, 197S]: 



with 



S4'k = J exp{ik • x)d^x . 

In the course of inflation, the interesting scales inflate outside the horizon and quantum fluctuations 
'freeze in' as classical fluctuations of the scalar field. This leads to energy density perturbations 
according to 

where V denotes the potential of the scalar field. By causality, on scales larger than the effective 
particle horizon, these perturbations cannot grow or decay by any physical mechanism. The spurious 
growth of super-horizon size density perturbations in certain gauges like, e.g., synchronous gauge, is 
a pure coordinate effect! 

This mechanism yields density perturbations Dx of a given size A which have constant amplitude 
at the time they re-enter the horizon after inflation is completed: 

Dx{t = X) = A. 

A natural consequence is thus the well known scale invariant Harrison-Zel'dovich spectrum 
iHarrison, 1970| , Zel'dovich, 1972| ]. Under certain conditions on, e.g., the potential of the scalar field 



(which may require fine tuning of coupling constants), one can achieve that these fluctuations have 
the required amplitude of A ^ 10^^. 

Recently, the failure to account for the largest scale structure of the otherwise so successful cold 
dark matter (CDM) model has caused some effort to find inflationary models with spectra that differ 
substantially from Harrison-Zel'dovich on large scales |Polarski and Starobinsky, 1992]. Although 



of principle interest, this work may turn out to be unimportant after the new COBE results have 
shown such a striking consistency with a scale invariant spectrum. 

B) Seeds as initial perturbations, topological defects: Initial fluctuations might have been 
triggered by seeds, i.e., by an inhomogeneously distributed matter component which contributes only 
a small fraction to the total energy density of the universe. Examples of seeds are a first generation 
of stars, primordial black holes, bosonic stars, cosmic strings, global monopoles and textures. We 
restrict our discussion to the latter three, the so called topological defects, which can arise naturally 
during phase transitions in the early universe [[Kibble, 1980| . 



To understand how they form, consider a symmetry group G which is broken by a scalar field (/) 
to a subgroup H at a temperature Tc. The vacuum manifold of the cooler phase is then generally 
given by Mq = G/H. Since the order parameter field cp (Higgs field) has a finite correlation length 
i {^h = horizon size) which is limited from above by the size of the horizon, the field varies in 
A^o if compared over distances larger than ^. If the topology of Mq is non-trivial, the scalar field 
can vary in such a way that there are points in spacetime where, by continuity reasons, (j) has to 



leave A^o and assume values of higher energy. This is the Kibble mechanism [Kibble, 19 



^Except for extended and hyper-extended inflation, the universe during the inflationary phase is a de Sitter universe, 
expansion is driven by a cosmological constant. 
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The set of points with higher energy forms a connected sub-manifold without boundary in four 
dimensional spacetime. The dimensionality, d, of this sub-manifold is determined by the order, r, of 
the corresponding homotopy group 7rj.(A^o): d = 3 — r, for r < 3. The points of higher energy of 
(p are often just called 'singularities' or 'defects'. 

For illustration, let us look at the simplest example, where A4o is not connected, 7ro(A4o) 7^ 0: 
At different positions in space with distances larger than the correlation length ^, the field can then 
assume values which belong to disconnected parts of TWq thus, by continuity, (j) has to leave 
A^o somewhere in between. The sub-manifold of points of higher energy is three dimensional (in 
spacetime) and called a domain wall. Domain walls are disastrous for cosmology except if they 
originate from late time phase transitions. 

A non-simply connected vacuum manifold, 7ri(A^o) 0) leads to the formation of two dimensional 
defects, cosmic strings. Domain walls and cosmic strings are either infinite or closed. 

If Mq admits topologically non-trivial (i.e. not shrinkable to a point) continuous maps from the 
two sphere, : ^ Mq , then 7r2(7Wo) 7^ 0, one dimensional defects, monopoles, form. 

Finally, continuous mappings from the three sphere determine 713. If 'it3{Mq) 7^ 0, zero dimen- 
sional textures appear which are events of non-zero potential energy. Using Derricks theorem 
^Derrick, 1964| ], one can argue that a scalar field configuration with non-trivial tt^ winding number 
(i.e. a texture knot) contracts and eventually unwinds, producing a point of higher energy for one 
instant of time. This type of defect is discussed in more detail in Chapters 4 and 5. 

Topological defects are very well known in solid state physics. Important examples are vortices 
in a super conductor or the vorticity lines in a super fluid. All four types of defects discussed above 
can also be found in liquid crystals see e.g. Chuang et al. [1991] and references therein. 

Depending on the nature of the broken symmetry, defects can either be local, if the symmetry 
is gauged, or global, from a global symmetry like, e.g., in the Peccei-Quinn mechanism. In the 
case of local defects, gradients in the scalar field can be compensated by the gauge field and the 
energy density of the defect is confined to the defect manifold which has a thickness given roughly 
by the inverse symmetry breaking scale. On the other hand, the energy density of global defects is 
dominated by gradient energy, with a typical scale given by the horizon size at defect formation. The 
extension of the induced energy density perturbation is thus about the horizon size at its formation. 
This leads to a Harrison-Zel'dovich initial spectrum. The difference to initial perturbations arising 
from inflation is that density fluctuations induced by topological defects are not Gaussian distributed. 

Local monopoles would dominate the energy density of the universe by far (Oa/ ~ 10^^ !!) and 
must be excluded. Local textures are not energetic enough to seed large scale structure. But global 
monopoles and global textures are quite promising candidates. Cosmic strings, both global or local, 
are interesting. One of the models which we introduce in the next section is based on local cosmic 
strings. 

We will find in Chapter 4, that the properties of the scalar field (/> other than the induced 
homotopy groups, e.g., the specific form of the Higgs potential, are of no importance for the defect 
dynamics. But the probability of defect formation via the Kibble mechanism might well depend on 
A^o- The symmetry breaking scale just determines the energy of the defects, i.e., the amplitude of 
the fluctuations. From simulations of large scale structure formation one finds that successful models 
requires T^G ~ m^/m^i ~ 10~^, which corresponds to a typical GUT scale of ~ lO^^GeV. This can 
also be understood analytically: From the perturbation equations derived in Chapter 2, we shall see 
that defect induced structure formation leads to initial perturbations with amplitudes A ~ IGvrGT^. 
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1.3.2 Linear perturbation theory 



Since the amplitudes of the initial fluctuations are tiny, at early times their evolution can be calculated 
within linear cosmological perturbation theory. 

Chapter 2 is devoted to a thorough discussion of this subject, we thus skip it here. The reader 
just interested in the results is referred to Chapter 3 where some important consequences from 
cosmological perturbation theory are discussed. 



1.3.3 N— body simulations 



On relatively small scales (up to about 20Mpc or so) perturbations become non-linear at late times. 
Structure formation must then be followed by non-linear numerical simulations, whose input is the 
perturbation spectrum resulting from linear perturbation theory. Since at the time when the non- 
linear regime takes over, the corresponding scales are much smaller than the horizon size, and since 
the gravitational fields and velocities are small, Newtonian gravity is sufficient for these simulations. 
For a realistic calculation of the process of galaxy formation, the hydrodynamical evolution of baryons 
has to be included. The different heating and cooling processes and the production of radiation which 
might partially provide the X-ray background and might reionize the intergalactic medium have to 
be accounted for. Furthermore, the onset of nuclear burning which produces the light in galaxies 
has to be modeled, in order to obtain a mass to light ratio or the bias parameter (which may well 
be scale dependent). 

The inclusion of the heating and cooling processes of baryons into N-body simulations is a very 
complex computational task and the results published so far are still preliminary [|Cen et al ., 1990|, 



Pen, 1992 , |Cen and Ostriker, 1992 , Cen et al., 1991]. Only very recently, the first calculations tak- 
ing into account nuclear burning have been carried out [pstriker, 1992 ] . High quality simulations 
which only contain collisionless particles have been performed to very good accuracy 



[ Centrella and Melott, 1982 , White, 1986 , Davis et al., 1985/87 |, but unfortunately they leave open 
the question how light traces mass and therefore, how these results are to be compared with obser- 
vations of galaxy clustering. 

I shall not discuss this important and difficult part of structure formation in this review and refer 
the reader to the references given above. 



1.3.4 Comparison with observations 

Of course the final goal of all the effort is to confront the results of a given model with observa- 
tions. The easiest and least uncertain part is to analyse the induced microwave background fluc- 
tuations. They can be calculated reliably within linear perturbation theory. Observations of the 
microwave background are also a measurement of the initial spectrum and it is very remarkable that 



the new COBE results are compatible with a Harrison-Zel'dovich spectrum [ Wright et al., 1992 
[Smoot et al., 1992| . 



On smaller scales {6 < 6°), the CMB anisotropics depend crucially on the question whether the 
universe has been reionized at early times {z > 100) and therefore indirectly on the formation of 
non-linearities in the matter density perturbation (which would provide the ionizing radiation). The 
different possibilities how the process of structure formation can induce anisotropics in the cosmic 
microwave background are discussed in detail in Chapter 3. 

On scales up to about 20Mpc, the up today most extensively used tool to compare models with 
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observations is the galaxy galaxy correlation function, ^gg('^)) 



1 



n{x)n{x + re)d xdVlg 



where n is the number density of galaxies. Usually, the amplitudes of initial fluctuations in a given 
scenario are normalized by J^{R = IQMpc): 

MR) = J ^GG{r)WR{rydr , 

with an observational value of Js{10Mpc)/{10Mpc)^ ~ 0.27. Here WR{r) is a (Gaussian or top hat) 
window function windowing scales smaller than R. Since the new COBE results now provide the 
amplitude of fluctuations on very large scales which are not influenced by non-linearities, future 
calculations will clearly be normalized on these scales. Observations fix the amplitude and slope of 
the galaxy galaxy and cluster cluster correlation functions to be 



iGc{r) 



r 



1.8 



\-1.8 



(1.5) 
(1.6) 



For galaxies one finds tq « 5 Ah ^Mpc [ Peebles, 1988 1, whereas for clusters ro depends on the 
"richness class" of the clusters considered. For rich clusters ro ~ (20 — 25)h~^ Mpc 



[ Bahcall and Soneira, 1983 1. Recently, a smaller amplitude for the cluster cluster correlation function 
has been found from the APM survey, S,cG = 4^GG and ^CG = '^CcG [ Dalten et al., 1991 ]. 

A disadvantage of the correlation function is its insensitivity to lower dimensional structures like 
sheets and filaments^. 

There are various other statistical tests which one can perform and compare with the sparse 
observations. I just mention a few: 

The Mach number, which gives the ratio between the average velocity and the velocity dispersion on 
a given length scale, M =< v >^ /o"^ 



Ostriker and Suto, 1990(| . 



The genus test, where one calculates the number of holes minus the number of islands in an iso- 



density contour | Gott et al., 1986 1. 

The 3-point or 4-point correlation functions, which determine the deviation from Gaussian statistics 
of the distribution of perturbations on a given length scale [ Peebles, 1980(] . They can be cast in the 
skewness, < 6p^ > / < ^P^ >^^^ and the kurtosis < 6p'^ > / < ^P^ 

Additional more qualitative results are: - The earliest galaxies must have formed at z > 5 but 
that there still must be substantial galaxy formation going on at z ~ 2 — 1. 

- From observations [ De Lapparent et al., 1986| , De Lapparent et al., 1988 1, 

[Geller and Huchra, 198£[ we can conclude that galaxies are arranged in sheetlike structures around 
seemingly empty voids. 

- Velocity observations have found large (~ 100/i~^Mpc) coherent velocity fields with 
{v) 500/cm/s. 

- A successful model should of course also obtain flat rotation curves of galaxies which have been 



observed with increasing accuracy since the seventies [Rubin, 1983]. 



^ Recently the amplitudes of the correlation function have been criticized to depend crucially on the sample 
size and thus to be physically meaningless. New analyses [Einasto et al., 198(:, Pietronero, 19871 [Dayis et al., 1988 



[Coleman and Pietronero, 1992 1 have shown that ro depends on the sample size, hinting that the distribution of galax 
ies may be fractal up to the largest scales presently accessible in volume limited samples, i?max ~ 30/i~^Mpc. If this 
objection is justified, the normalization procedure with the help of the correlation function is useless! 
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1.4 Models 



The most simple picture, a universe with Qh'^ = fi^/i^ ~ 0.1 and adiabatic initial fluctuations 
is definitely ruled out by the limits of the microwave background fluctuations [ |Gouda et al., 1989 



pouda et al., 1991| ]. This is a first, very important result in the discussion of different models which 



might account for the formation of large scale structure. 

To give the reader a taste of the presently favored scenarios, we present here five cases. A sixth 
possibility, the "texture scenario", will be discussed in Chapters 4 and 5. It is clear that a mixture of 
the hot dark matter (HDM) and cold dark matter (CDM) models presented here, as well as defects 
with HDM or any of the proposed scenarios with addition of a cosmological constant might lead to 
models that fit the observations better. The ones presented here are partly chosen by reasons of 
simplicity. The first attempt to treat a class of models systematically is given by Holtzman [1989], 
where 94 different combinations of {Q\,^Icdm,^hdm,^b) are investigated. 

In this section we give a short description of the models chosen, and compare them with some 
observational phenomena. Our aim is only to give a sketchy overview of these models; readers 
interested in more details are referred to the literature given in the text. 

1.4.1 The isocurvature baryon model 

As the most conservative alternative to the adiabatic baryon model, Peebles [1987] pursued the 
question, whether it is possible to construct a viable scenario of structure formation without the 
assumption of any exotic, i.e. up to now unobserved, form of energy, a universe with J7o = ~ 0.1, 
which is still marginally consistent with nucleosynthesis limits on 0^. 

In order not to overproduce CMB fluctuations, one has to assume isocurvature initial perturba- 
tions, i.e. no perturbations in the geometry on cosmologically relevant scales (see Section 3.1). From 
this one can conclude 6p/p ~ on scales which are larger than the size of the horizon. Since the 
universe is radiation dominated initially, pr ^ Pb this yields \5pr/pr\ \Spb/pb\, i-e. isothermal 
fluctuations. 

Isocurvature perturbations allow for relatively high initial values of 5pB and therefore lead to early 
structure formation. Galaxies form at Zg > 10. The even earlier formation of small objects reionizes 
the universe. Small scale anisotropics in the CMB can then be damped by photon diffusion. For 
photon diffusion to be effective, reionization must take place before the universe becomes optically 
thin. This provides a lower limit to a 'useful' reionization redshift (see Chapter 3): 

^i>^dec~ 100(^)2/3 . 

Photon diffusion then damps fluctuations on all scales smaller than the horizon scale at decoupling, 
In^Zdec) and correspondingly all angular scales with 



e < Odec ~ 6y^^(100/zrfec) . 

The problems of this scenario are twofold: First, the quadrupole anisotropy of the CMB turns out to 



be unacceptably large [ Gouda et al., 1991 |. A way out of this problem is a steep initial perturbation 
spectrum, but then it is difficult to reproduce the large amplitudes of the galaxy correlations on 
large scales. (Although, for steep a enough spectrum, even quadrupole fiuctuations may be damped 
by photon diffusion, see Section 3.2.3!) 

Secondly, observations hint that galaxy formation might peak around z ~ 2 or, at least, is still 
going on around z ~ 1, whereas in the isocurvature baryon model the process of galaxy formation is 
most probably over before 2 = 3. 
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1.4.2 Hot dark matter 



Massive neutrinos with "^i ~ 200eVh^ they can provide the dark matter of the universe and 



dominate the total energy density with 1 = ^Itot ^ [ Doroshkevich et al., 1980, Bond et al., 1980 1. 
Large scale structure then develops as follows: 

Initial fluctuations from inflation give rise to a scale invariant spectrum of Gaussian fluctuations. 
These initial fluctuations are constant until they 'reenter the horizon' (i.e. their scale becomes smaller 
than the size of the horizon, I < Ih). Thereafter they decay by free streaming if the universe is still 
radiation dominated; if the universe is matter dominated {z < 2 x 10^) they grow in proportion to 



the scale factor [ Bond and Szalay, 1983 , [Durrer, 1989a[| . This leads to a short wavelength cutoff of 



the linear perturbation spectrum at Ips ~ 'iO{m,y /30eV)~^ Mpc. The corresponding cutoff mass 
is Mps ~ W^^{30eV/mu)'^MQ. The linear spectrum for HDM is given in Fig. 3. 

In this model, large objects with mass ~ Mps (large clusters) form first . They then frag- 
ment into galaxies. Gravitational interaction of collisionless particles generates sheets, pancakes 
[iZel'dovich, 197C ], and galaxies are thought to lie on the intersections of these sheets. This leads to 



a filamentary structure. Simulations of HDM show [Centrella and Melott, 1982, White et al., 1983 



that in order to obtain the correct galaxy correlation function today, the large scale structure be- 
comes heavily overdeveloped (see Fig. 4). The other main problem is that galaxy formation starts 
only very recently (z ~ 1). The model has serious difficulties to account for quasars with redshifts 
z > 3 — 4. In addition to these grave objections, CMB fluctuations turn out to be too large in this 
model. Since galaxy formation is only a secondary process, initial fluctuations which determine the 
amplitudes of CMB anisotropies must be rather large. 

Because their thermal velocities are relatively high, massive neutrinos can only marginally provide 
the dark matter of galaxies but they cannot be bound to dwarf galaxies. If neutrinos are to constitute 
the dark matter of a virialized object with velocity dispersion a and size r, their mass is limited from 



below by the requirement |Tremaine and Gunn, 1979 1 



n.. > 30ey(^^^)V^(-^)V^ . 

a Wkpc 

Since also dwarf galaxies do contain substantial amounts of dark matter 



[Carignan and Freeman, 1985], this is another serious constraint for the HDM model. 



1.4.3 Cold dark matter 

In this scenario one assumes the existence of a cold dark matter particle which at present dominates 
the universe with 1 = ^Itot ~ ^cdm- Particle physics candidates for such a matter component are 
the axion or the lightest super-symmetric particle. A more extended list can be found in Kolb and 
Turner [1990]. 

Again it is assumed that an inflationary phase leads to a scale invariant spectrum of Gaussian 
initial fluctuations. After the universe becomes matter dominated at Z(.q ~ 2 x 10^ /i^, perturbations 
smaller than the horizon, I < leg ^ 10(J7/i2)-i. start growing. Damping due to free streaming is 
negligible. Because of the logarithmic growth of matter fluctuations in the radiation dominated era 
(see Section 3.1), the spectrum is slightly enhanced on smaller scales / < /eg- 

According to this linear spectrum, sub-galactic objects form first. Once the perturbations become 
non-linear, these objects virialize and develop flat rotation curves. Only very recently, big structures 
begin to form through tidal interactions and mergers pavis et al, 1988|] . 
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To obtain a mass distribution with i^dyn ^ 0.1 — 0.2 < i7 = 1, it is necessary that most of the 
mass is in the form of a dark background which is substantially less clustered than the luminous 
matter. This can be achieved with the idea of biasing: Luminous galaxies only form at high peaks 
of the density distribution. Since for a Gaussian distribution high peaks are more strongly clustered 



than average, this simple prescription has the desired effect [ Kaiser, 1985 , Bardeen et al., 1986 1. 
Usually one introduces a bias parameter b and requires that galaxies form only in peaks of height 
6(7, where a is the variance of the Gaussian distribution of density peaks. The best results are 
obtained for a bias parameter 6 ~ 1.5 — 2. Clearly, once the correct hydrodynamical treatment 
of baryons is included in the numerical simulations, such a bias parameter (which otherwise is just 
an assumption of how light traces mass) could be calculated. First preliminary results of such 



calculations [ Cen et al., 1990 , Ostriker, 1992 ] indicate that the above assumption is probably quite 
reasonable. A problem of the biasing hypothesis is the prediction of many dwarf galaxies in the voids 
which have not been found despite extensive searches. 

The large scale structure obtained in this scenario looks at first sight rather realistic (see Fig. 5). 
It leads to the right galaxy galaxy correlation function up to lOMpc. But the recently detected huge 
structures like the great wall (see Geller and Huchra [1989]) are very unlikely in this model. This 
situation has been quantified by comparing the angular correlation function from the deep IRAS 



survey with the one predicted by CDM [ Maddox et al., 1990 [. There, a substantial excess of power 
(as compared to the CDM model) on large angular scales, 9 > 2° is found (see Fig. 6). 

Since galaxies form relatively late (at z ~ 2), it might also be difficult to produce the very high 
redshift quasars (z ~ 5). But since their statistics are still so low, and since very little is known on 
the ratio of normal galaxies to quasars at high redshift, this may not be a real problem. 



1.4.4 Cosmic strings 

Here, initial fluctuations are seeded by cosmic strings which form via the Kibble mechanism (see 
Section 3) after a phase transition at Tc ~ W^^GeV [Vilenkin, 198C[. 

Inter-commutation and gravitational radiation of cosmic string loops (see Vachaspati and Vilenkin 
[1984], Durrer [1989b]) determine the evolution of a network of (non-superconducting) cosmic strings. 
Numerical simulations support analytical arguments for a scaling law, Pstring oc l/{at)'^ oc p, 
for the energy density of a cosmic string network. In contrast to gauge monopoles or domain 
walls, strings do not dominate the energy density of the universe and are cosmologically allowed 



[lAlbrecht and Turok, 198g| ]. 

It is well known that a static straight cosmic string does not accrete matter, whereas a cosmic 
string loop from far away acts like a point mass [ Vilenkin, 1980 [. High resolution simulations of 
cosmic string networks ]Bennett and Bouchet, 1990 [ have shown that the loops that chop off the 
network are too small {hoop < 0.01///) for efficient accretion. But moving long strings produce large 
accretion wakes behind them which might provide the sheets and walls observed in the universe 



Bertschinger, 1987 , Perivolaropoulos et al., 199C ]. In this scenario galaxies form via fragmentation 



and/or accretion onto loops. 

A relatively new idea is that chopping off small loops could lead to wiggles on the long strings. 
Cosmic strings with such small scale wiggles give rise to strings with an effective tension which is 
smaller than its effective energy density. In contrast to the original cosmic strings, wiggly strings can 
accrete matter even if they are static [ parter, 1990 , Vilenkin, 199C, [Vachaspati and Vilenkin, 1991 ]. 
In order for the large scale structure not to be 'drowned', in these small scale structures, the scenario 
works best if the dark matter is hot, massive neutrinos (HDM), so that the small scale fluctuations 
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are damped by free streaming. 

The requirement for successful structure formation on one hand and the difficulty not to over- 



produce microwave background anisotropies [Btebbins, 1988| on the other hand tightly constrain the 
possible value of the symmetry breaking scale 

^ = G7?2 « (1 - 2) X 10-^ . 

Recent work using new cosmic string simulations in which these estimates were redone to compare 
the CMB anisotropies with the COBE measurements led to a similar value 
[Bennett and Bouchet, 1992, Perivolaropoulos, 1993]. 



1.4.5 Global monopoles 

Like for cosmic strings also the energy density of global monopoles produced by the Kibble mech- 
anism obeys a scaling law. Therefore, they are candidates for a model of structure formation. In 
contrast to local monopoles, the gradient energy of global monopoles introduces a long range in- 
teraction, so that monopole anti-monopole pairs annihilate, leaving always only a few per horizon 



[Barriola and Vilenkin, 1989, Bennett and Rhie, 1990] . 

One assumes, like in the CDM model, that the matter content of the universe is dominated by cold 
dark matter with 1 = 17 QcDM- The large scale structure induced by global monopoles seems to 
look quite similar to the texture scenario. Galaxy formation starts relatively early, the galaxy correla- 
tion function and large scale velocity field are in agreement with observations. The CMB fluctuations 
are similar to those obtained for the texture scenario (see Chapter 5) [ [Bennett and Rhie, 1992 1. 



Recent investigations [ Pen et al., 1993| ] claim that all models with global defects and CDM miss 
some large scale power on scales / > 20Mpc. 
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Chapter 2 

Gauge Invariant Perturbation Theory 



For linear cosmological perturbation theory to apply, we must assume that the spacetime manifold M. 
with metric g and energy momentum tensor T of "the real universe" is somehow close to a Friedmann 
universe, i.e., the manifold Ai with a Robertson- Walker metric g and a homogeneous and isotropic 
energy momentum tensor T. It is an interesting, non-trivial unsolved problem how to construct g 
and T from the physical fields g and T in practice. There are two main difficulties: Spatial averaging 
procedures depend on the choice of a hyper-surface and do not commute with derivatives, so that 
the averaged fields g and T will in general not satisfy Einstein's equations. Furthermore, averaging 
is in practice impossible over super-horizon scales. 

We now assume that there exists an averaging procedure which leads to a Friedmann universe 
with spatially averaged tensor fields Q, such that the deviations (T^;/ — T^,^)/ max^afi} {\Tai3\} and 
idfiiy ~gfiiy)/ ''^(^^{affjl'Sas} stTe Small, and g and T satisfy Fricdmann's equations. Let us call such an 
averaging procedure 'admissible'. There might also be another admissible averaging procedure (e.g. 
over a different hyper-surface) leading to a slightly different Friedmann background (^2, ^2). In this 
case, the averaging procedures are isomorphic via an isomorphism (f) on M which is close to unity: 

92 = (t>*9 

where 4>*{Q) denotes the pushforward of the tensor field Q under cp. The isomorphism (p can be 
represented as the infinitesimal flow of a vector field X, cj) = (j)^ . Remember the definition of 
the flow: For the integral curve 7a; (s) of X with starting point x, i.e., jxis = 0) = x we have 
(t>^{x) = 7x(s). In terms of the vector field X, the relation of the two averaging procedures is given 



In the context of cosmological perturbation theory, the isomorphism is called a gauge transforma- 
tion. And the choice of a background (g, T) corresponds to a choice of gauge. The above relation is 
of course true for all averaged tensor fields Q and Q2- Separating Q into a background component 



by 



92 = <i>*g = g- eLxg + O(e^) 



(2.1) 
(2.2) 




(2.3) 
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Since each vector field X generates a gauge transformation (f> = cf)^ , we can conclude that only 
perturbations of tensor fields with LxQ = for all vector fields X, i.e., with vanishing (or constant) 
'background contribution' are gauge invariant. This simple result is sometimes referred to as the 
'Stewart Lemma' [Stewart and Walker, 1977]. 

The gauge dependence of perturbations has caused many controversies in the literature, since 
it is often difficult to extract the physical meaning of gauge dependent perturbations, especially on 
super-horizon scales. This has led to the development of gauge invariant perturbation theory which 
we are going to present in this chapter. The advantage of the gauge-invariant fromalism is that the 
variables used have simple geometric and physical meanings and are not plagued by gauge modes. 
Although the derivation requires somewhat more work, the final system of perturbation variables is 
usually simple and well suited for numerical treatment. We shall also see, that on subhorizon scales, 
the gauge invariant matter perturbations variables approach the usual, gauge dependent ones, and 
one of the geometrical variables approaches the Newtonian potential, so that the Newtonian limit 
can very easily be performed. 

There are two review articles on the subject |Kodama and Sasaki, 1984| , [Mukhanov et al., 1992| |. 
Our treatment will be in some way complementary. Collisionless particles, photon propagation and 
seeds (Sections 2.3, 2.4, 2.5) are not discussed in these reviews. On the other hand, we do not 
investigate perturbation theory of scalar fields which is presented extensively in the other publica- 
tions, since it is needed to treat fluctuations induced by inflation. We want to discuss perturbations 
induced by topological defects. Here, the scalar field itself is a small perturbation on a matter or 
radiation dominated background. This issue will become more clear later. 

First we note that since all relativistic equations are covariant (i.e. can be written in the form 
(5 = for some tensor field Q), it is always possible to express the corresponding perturbation 
equations in terms of gauge invariant variables. 



2.1 Gauge Invariant Perturbation Variables 

In this section we introduce gauge invariant variables which describe the perturbations of the metric 
and the energy momentum tensor in a Friedmann background. 

There are two main approaches to find gauge invariant quantities: One possibility is to make 
full use of the above statement that tensor fields with vanishing background contribution are gauge 
invariant and use them to define gauge invariant perturbation variables. Examples are the Weyl 
tensor, the acceleration of the energy velocity field, anisotropic stresses, shear and vorticity of the 
energy velocity field. This covariant approach was originally proposed by Hawking [1966] and later 
extended by Elhs and Bruni [1989], Elhs et al. [1989], Hwang and Vishniac [1989], Bruni et al. 
[1992b], Dunsby [1992] and others. 

The other possibility is to arbitrarily parametrize the perturbations of the metric, the energy 
momentum tensor, the distribution function, a scalar field... and discuss the transformation proper- 
ties of these gauge dependent variables under gauge transformations. One can then combine them 
into gauge invariant quantities. This way was initiated by Gerlach and Sengupta [1978] and Bardeen 
[1980] and later continued by Kodama and Sasaki [1984],Kasai and Tomita [1986], Durrer [1988], 
Durrer and Straumann [1988], Durrer [1989a], Durrer [1990], Mukhanov et al. [1992] and others. 
In this approach one usually performs a harmonic analysis and the gauge invariant perturbation 
variables found in this way may be acausal (e.g., they may require inverse Laplacians over spatial 
hyper-surfaces) . 

As in the second approach we divide the perturbations into scalar, vector and tensor contributions. 
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but we do not perform the harmonic analysis. Our perturbation variables are thus space (and time) 
dependent functions and not just amplitudes of harmonics. In order to obtain unique solutions, we 
require (for non-positive spatial curvature) all the perturbation variables to vanish at infinity. I 
partly relate the two approaches by originally performing the second, but in many cases identify 
gauge invariant variables with tensor fields which vanish in an unperturbed Friedmann universe. It 
could not be avoided to use a somewhat extensive vocubulary for all the variables used in this text. 
To help the reader not to get lost, I have included a glossary in Appendix B. 
Like in Chapter 1, the unperturbed line element is given by 

ds^ = a'^{—dt^ + ^ijdx'^dx^) , 

where 7 is the metric of a three space with constant curvature A; = ±1 or and overdot denotes 
derivatives w.r.t conformal time t. We first define scalar perturbations of the lapse function a, the 
shift vector /3 and the 3-metric g of the slices of constant time|^ by 

a = a{l + A) (2.4) 

(3 = aB^'di (2.5) 

g = a^[{l + 2HL-{2/3)l'^AHThij + 2l'^HT\ij]dx'dx^ . (2.6) 

We denote 3-dimensional vector and tensor fields by bold face letters; | and A denote the covariant 
derivative and Laplacian with respect to the metric 7. The variables A, B, Hl, and Ht are arbitrary 
functions of space and time. To keep them dimensionless, we have introduced a length / which, in 
applications, will be chosen to be the typical scale of the perturbations so that, e.g., 0{B) ~ 0{lB^^). 

By choosing the metric perturbations in the form (2.4-p.6|), we restrict ourselves to scalar type 
perturbations, but we do not perform a harmonic analysis. Vector perturbations of the geometry are 
of the form 

f3 = aB'di (2.7) 

g = a'^lj-j + lHiy + lHj^i]dx'dx^ , (2.8) 

where B^di and H^di are divergencefree vector fields which vanish at infinity. 
Tensor perturbations are given by 

g = a'^[-y.. + 2Hij]dx'dx^ . (2.9) 

Here Hij is a traceless, divergencefree, symmetric tensor field on the slices of constant time. 
Writing the 4-dimensional metric in the form 

gfiu = g^iv + a^h^y , (2.10) 

the above definitions of the perturbation variables for scalar perturbations yield 

hiS) = -2A{dtf + 2lB,i dtdx' + 2{Hl AHrhijdx'dx^ + 2f HT\,jdx' dx^ . (2.11) 

3 

For vector and tensor perturbations one obtains correspondingly 

/i(^) = 2Bidx'dt + l{Hi\j + Hj\i)dx'dx^ and (2.12) 

h^^^ = 2Hijdx'dx^ . {2 A3) 

short description of the 3+1 formahsm of general relativity is given in Appendix A 
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From (2.4), ( p.5D , ( |2.6D , (2.7) and (2.8), one can calculate the 3-dimensional Riemann scalar and 
the extrinsic curvature given by Kij = —rii-j. Palatini's identity (see Straumann [1985]) yields the 
general formula 

6R{g + Sg) = {6gl)l^-A6gl. 

This leads to the following result for scalar perturbations 

72 

6R= -4:a-'^{A + 3k)n , 71 = Hl- -AHt , (2.14) 

j^ianiso) ^ _^2al^(^HT\ij - ^^ijAHr) + a{laiij - ^--fij Aa)]dx' dx^ , la = I^Ht - IB . (2.15) 

For vector perturbations one obtains 

6R^ = (2.16) 

j^V{amso) ^ _^a{Hi\j + Hj\i) + '^{ai\j + aj\i)\dx'dx^ with ai = IHi - Bi . (2.17) 



j^{amso) jg traceless contribution to the extrinsic curvature of the slices of constant time or, what 
amounts to the same thing, the shear of the normal to the slices. 

We now investigate the gauge transformation properties of the variables defined above. We 
introduce the vector field describing the gauge transformation by 

X = Tdt + Vdi . 

Using simple identities like Lx{df) = d{Lxf) and {Lxj)ij = ^i\j+^j\'h we obtain the Lie derivative 
of the unperturbed metric 

Lxg = a^[-2(-T + f)dt^ + 2(Li - Tu)dx'dt + (2-T-/ii + L,,,- + L.u)dx'dx^ . 
a ' a 

For scalar type gauge transformations, L* is of the form L* = /LI* for some scalar function L. Inserting 
this above and comparing with the parametrization of the perturbed metric h^^\ eq. ( |2.11| ), we find 
the following transformation laws: 

A A + -T + t (2.18) 
a 

IB IB-T + IL (2.19) 

Hl ^ Hl + -T + {1/3)AL (2.20) 
a 

I^Ht fHr + IL . (2.21) 

(2.22) 

This yields the transformation properties for IZ and a 

TZ^n+^T; a^a + T. (2.23) 
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The following scalar perturbation variables, the so called Bardeen potentials, are thus gauge invariant 
(see [ Bardeen, 198C1| , Kodama and Sasaki, 1984 1 and [ Durrer and Straumann, 1988| ): 

$ = 7^ - (a/o)/f7 (2.24) 
= A- {a/a)la - la . (2.25) 

We shall later see that is a relativistic analog of the Newtonian potential. 

For vector type gauge transformations, where T = and L* is a divergence free vector field on 
the hyper-surfaces of constant time, one immediately sees that ai is gauge- invariant. 

Tensor perturbations are of course always gauge invariant (there are no tensor type gauge trans- 
formations) . 

Gauge transformations remove two scalar and one vector degrees of freedom, so that geometrical 
perturbations are fully characterized by six degrees of freedom which we parametrize in the gauge 
invariant variables ^,^,ai and Hij. 

Some combinations of these quantities have simple geometric interpretations. From the Stewart 
Lemma, we know that the Weyl tensor, which vanishes in an unperturbed Friedmann universe, is 
gauge invariant. A somewhat lengthy calculation shows that the electric and magnetic components 
of the Weyl tensor, are given by 



Eij — —Ciojo 



1 

a 

-[{<!>- - - 



[(CD - ^.)|,^. - -A(cl> - x&)7,, + a,,] (2.26) 



= <^{ilrrMj)m\l " ■^lj)l(^m\s\ {2.27) 

where = + and (i..j) denotes symmetrization in i and j 

is the sum of vector and tensor contributions to the extrinsic curvature (this result is obtained in 
Bruni et al. [1992b]). 

We now discuss perturbations of the energy momentum tensor. We define the perturbed energy 
density p and energy velocity field u as the timelike eigenvalue and eigenvector of the energy mo- 
mentum tensor (note that, apart from symmetry, we do not make any assumptions on the nature of 

^ w = —pu , U = —1 . 
We then define the perturbations in the energy density and energy velocity field by 

p = p{l + 6) , (2.28) 
u = u^dt + u'di ; (2.29) 

is fixed by the normalization condition, = a~^(l — A). In the 3~space orthogonal to u we define 
the stress tensor by 



'^'^ = P^^P^T'^P , (2.30) 
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where P = u(i$ u + g is the projection onto the sub-space of TqM. normal to u{q). One obtains 

Tq — Tj — Tq — u . 

The perturbations of pressure and anisotropic stresses can be parametrized by 

t/ =p[{1 + ttl)6/ +7ri^] ,with ^1 = 0. (2.31) 

Again, we decompose these perturbations into scalar, vector and tensor contributions. For scalar 
perturbations one can set 



u 



- (1 - A) , - = -Iv'^ and vr^ = I^tt^]^ - -Avr^^^ 



As before, the Lie derivatives of the unperturbed quantities p , u = a ^dt , r = pdi dx^ 
determine the transformation laws of the perturbation variables. One finds 

Lxp = Tp = -3{l + w)^Tp (2.32) 

Lxu = [X, u] = -a-i [(-T + t)dt - L'di] (2.33) 

a 

hxf = Lxpdi (g) dx' = Tpdi ® dx* - 3— (1 + w)-Tpdi ® dx' , (2.34) 

w a 

where we have used the background energy equation, 

- \ ^ - 

p = -3(1 + w)-p , 
a 

and the definitions 

w = p/p , cl= p/p . 
From eqs. (|2.32D - (|2.34| ) we obtain the transformation laws 

6 ^ 6- 3(1 + w)-T 
a 

Iv Iv + IL 

2 

VTL ^ TTL -3^(1 + ?i;)-r 
w a 

vr — > vr . 

A first gauge invariant variable is therefore 
n = vr , 

the scalar potential for anisotropic stresses. The other gauge invariant combination which can be 
constructed from matter variables alone is 

T = TTL- {cl/w)5 . 
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Defining an entropy flux 5^^ of the perturbations in the sense of small deviations from thermal 
equilibrium [ Straumann, 1985 , Appendix B ], one finds for the entropy production rate of the 
perturbation [Durrer and Straumann, 1988 1 

5'^;^=3(a/a2r)r , 

where T denotes the temperature of the system. The variable T thus measures the entropy production 
rate. 

We now split the covariant derivative of the velocity field in the usual way into acceleration, 
expansion, shear and vorticity: 

Here 

a = VuU is the acceleration, 
denotes the projection onto the sub-space of tangent space normal to u, 6 = u'^^ is the expansion, 
is the shear of the vector field u and 



1 



CO 



lP^PP[uX;p-Up,x) 



2 



is the vorticity. A short calculation shows lo^i, = 0, cro^t = and 

= jViij with IV = lv- PHt ■ (2.35) 

We choose V as gauge invariant scalar velocity variable. For the acceleration one obtains = and 

ai = A,i = +Vi + -Vi , 
a 

which shows for the first time the analogy of ^ with the Newtonian potential. 

There are several different useful choices of gauge invariant density perturbation variables: 

Ds = 5 + 2,{l + w){a/a)la (2.36) 

Dg = 5 + 2,{l + w)n = Ds + 'i{l + w)^ (2.37) 

D = Ds + 2,{l + w){a/a)lV . (2.38) 

For a physical interpretation of these variables note that 

= Ptp,f. (2.39) 
(!),)(,,•)+ 3(1 + «;)M/(,,) = PtP^P■,^.. (2.40) 
(Z?,)(,,)+ 3(1 + u;)(vl/ -<!>)(,,) = PtPj'P;f.u. (2.41) 
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Here we have set S^^jj = S^ij — (l/3)A5'7ij for an arbitrary scalar quantity S. 

Therefore, D and Dg + 3(1 + w)'^ are potentials for the first and second "spatial derivatives" of the 
energy density. 

For vector perturbations only 

u' = -v' and7r^. = ^(4 + ^,'') 

survive. Vector type gauge transformations yield the transformation laws 
v' v' + L' 

it' tt' . 

In addition to the anisotropic stress potential tt^ = W , two interesting gauge invariant quantities 
are the shear and vorticity of the vector field u: 

Ui-j + Uj-i = a{V,\j + Vj\i) , with V' = v' - IW , (2.42) 

Ui-j - Uj-i = a{VLi\j - Q.j\i) , with = - . (2.43) 

For tensor perturbations the only variable tt^j = Wj is of course gauge invariant. 

We now show that for perturbations which are small compared to the horizon distance, /j^, in a 
generic gauge the gauge invariant combinations V and approach the original v and S. Let us 
choose our free length scale I to be the typical size of a given perturbation. From the above equation 
it is then clear that for / ^ /// = t, D ^ Dg. 

Noting that perturbations of the Einstein tensor are given by second derivatives of the metric 
perturbations (Palatini's identity, see e.g. Straumann [1985]), we obtain the following order of 
magnitude equation: 

0(—)0{87rGT.,) = 0{t-^^-^ + (/i)"^^ + T^^) . (2.44) 
T 9 9 9 

Using Friedmann's equation 

0(87rGT^,) = 0{a/af = 0(1/*^) = 0{l/l]j) 
this yields 

o(f) = o(^ + (WO^ + (WO^^). 

T 9 9 9 

On sub-horizon scales, Ih 3> I the metric perturbations are thus generically much smaller than the 
matter perturbations and the difference between the gauge invariant quantities V, -D^ ), V^, \Orrf 
and V, (5, v'^ becomes negligible. 



2.2 The Basic Equations 

In this section we write down the perturbation equations resulting from Einstein's equation, and 
energy momentum "conservation" in a form which will be convenient later. All these equations are 
most easily derived using the 3+1 formalism of gravity (see Appendix A) as we shall demonstrate 
for a few examples. 
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The perturbations of Einstein's equations and energy momentum conservation can be expressed 
in terms of the gauge invariant variables defined above. (A simple derivation of the equations for 
scalar perturbations is given in Durrer and Straumann [1988].) To simplify the notation we now 
drop the bar over background density and pressure. 

A) Constraint equations 

47rGaV-C = -(A + 3A;)$ (2.46) 
4.T:Ga^{p + p)lV = (d/a)^' - 6 . (2.47) 



B) Dynamical equations 

-87rGaV^An = A($ + *) (2.48) 
87rGa2p(r + {cl/w)Dg - (2/3)^^ AH) = (d/a){^ - [(l/a)(^)T} + 

{2a(d/a2)- + 3(d/a2)2}[*- l/a(— )■] . (2.49) 

a 

Since vector and tensor type perturbations are not treated in Durrer and Straumann [1988] and 
Durrer [1990], we present an explicit derivation of the vector perturbation equations, making use of 
the 3+1 formalism of general relativity (see Appendix A and Durrer and Straumann [1988]). For 
vector perturbations, the unit normal to the equal time slices is given by 

n = a-\dt-l3) = a-\dt-B'di) . 

We now decompose the energy momentum tensor in the form 

T = €n(^n + n<^S + S(^n + T 

where, as before, bold type vector and tensor fields, S and T, are tangent to the equal time hyper- 
surfaces. Using the Gauss-Codazzi-Mainardi formulas to express the four dimensional curvature in 
terms of the three dimensional and the second fundamental form, one can derive the following 3+1 
split of Einsteins equations (Appendix A4) : 

R+{trKf -triK"^) = WnGe (2.50) 

V-K-Vtr{K) = SttGS (2.51) 
dtK - LpK + Hess(a)- 

-a[mcci - 2K ■ K + {trK)K] = -ATrGa[2T + g{e - trT)] , (2.52) 

where K = -^{L^g — dtg) is the second fundamental form. 

For vector perturbations only the second constraint equation and the traceless part of the dy- 
namical equation contribute. From our definiton of vector perturbations of the energy momentum 

tensor, one finds 

S = a{p + p){v'-B')di . 
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In what follows we use the notation X^J^ for the symmetric, traceless tensor constructed from the 
vector field X, i.e., X^J^ = ^{X^j + Xj\i). We then have 

g = a'^{-fij + 2lH^p)dx'dx^ 
and with ( |2.17| ) the second fundamental form is given by 



This leads to 
V K = 



1 

2^ 



1 

2^ 



;A + 2k)(7'di 



The constraint equation ( 2.51 ) thus results in 
1 



-(A + 2k)a' = 8ttG{p + p)a^Ve 



(2.53) 



Let us now proceed to the dynamical equation. Up to first order we obtain the following expres- 
sions for the terms in (2.52): 



Ricci{g) 

SttGT 
47rGgi(e - trT) 
dtK 



-2dB^pdx'dx^ 

''J 

[(«)2( + 2///(f )) + 2^alJWdx^ 

2k{-fij + 2lHip)dx'dx^ 

SnGa^pi-fij + 21h\P + lT:f^^)dx'dx^ 

AttG{p - 3p)a2(7ij + 2lH^p)dx'dx^ 

-[a(7., + 21hI^^) + d{a^^ + 2lHiP) + aaiP]dx'dx^ 



The result for Ricci{g) is again easily derived using Palatini's identity [ Straumann, 1985 ]. 
With the help of the background relation 

A7rGa^{p-p) = - + i-)^ + 2k , 
a a 

equation ( |2.52 ) then yields 

+ 2(^).f =8vrGaVnSr^ 

Since we require lim^^oo = lim^^oo = 0, the tensor fields, a^J^ and ir^J^ uniquely determine 
the corresponding vector fields: 



ai + 2(-)cjj 
a 

or (a^cr,:)- 



S-rrGa'^pmi 
SirGa^plUi 



(2.54) 
(2.55) 



In the absence of anisotropic stresses, vector anisotropics in the extrinsic curvature thus decay like 
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In a similar way one finds the tensor perturbation equation 

iiij + '^^Hij + {2k - A)Hij = SfrGa'^pUij . (2.56) 

This is a wave equation with source term. It describes the creation, propagation and damping of 
gravitational waves in a Friedmann background . 

C) Conservation equations 

The energy and the momentum conservation equations of each independent type of matter (i.e. 
each matter component which does not interact other than gravitationally with the rest) yields the 
following equations of motion for the scalar perturbation variables Da and Va, where the index „ 
denotes the different matter components (e.g. radiation, dark matter, baryons ...): 

Da-3wa{d/a)Da = {A + 3k)[{l + Wa)lVa + 2{a/a)wJ^Ua] 

+3{l + Wa)47rGa\p + p){lV -IVa) , (2.57) 

2 

lVa + {d/a)lVa = —^Da + -^^ra + '^ + 2/3{A + 3k)-^^fUa. (2.58) 
The total perturbations are defined as the sums: 

pD = Y,PaDa , ip+p)V = J2iPc.+Pa)Va . (2.59) 

a a 

The adiabatic sound speed, Ca and enthalpy, Wa are 

Ca=Pa/Pa , Wa=Pa/Pa and W = = p/p , cl=p/p. 

l^a Pa 

The corresponding equations for interacting matter components are derived in 



[Kodama and Sasaki, 1984]. 



We shall later also use equations (2.57) and (2.58) for a one component fluid in terms of the 
density perturbation variable Dg and V. One easily finds 

Dg + 3(c2 - w){d/a)Dg - (1 + w)lAV + 3?i;(d/a)r = (2.60) 

IV + (d/a)(l - 3c^)/y = (^ - 3c2$) + -^Dg + -^[T + |(A + 3k)fU] (2.61) 

1 + w 1 + w 6 

Again we derive the conservation equation for for vector perturbations in some detail. We start 
with the spatial part of the 3+1 split of 
T^^^ = (see Appendix A.3): 

-{dt - Lb)S = -V(lna)e + 2KS + {trK)S - -V • (aT) . (2.62) 

a a 

All the terms in this equation are readily calculated. One obtains for each matter component 
(suppressing the index a) 

S = {p + p)a-\v' -B')di 

-V(aT) = V • T = -^p(Am' + 2km')di 
a 2a^ 
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2KS = -2^{p + p){v' -B')di 

(trK)5 = -^^{p + p){v' -B')di 

Lf^S = (in first order) 

dtS = [-ha-\p + p){v' -B')+a-\p + p){v' -B')+a-\p + p){e -B'W 

= a-\p + p)[v,-Bi-{^ + 2,cl)'^{vi-Bi)]di . 

For the last equality sign, we made use of tlie background equation p + p = — 3(d/a)(l + c^)(/o + p). 
Inserting all these results, eq. ( |2.62| ) becomes 

tr + {I - 2,cl)-Q.' = -^—{l\ + 2k)m'. (2.63) 

a 2{p + p) 

If there are no sources present, 11* = , and if = to = constant, the amplitude of the vorticity is 
proportional to a^'^'>~^ . In comparison to the expansion velocity, d/a, the vorticity behaves like 

^^o.5(9«,-i) (2.64) 

a 

(as long as curvature is negligible, i.e. for k <C 1/i^)- Especially, an initial vorticity in a radiation 
dominated universe {w = 1/3) grows relative to the expansion velocity in the course of expansion. 

In addition to Einstein's field equation and the conservation equations which are of course a 
consequence of them, we have to add matter equations to fully describe the system. If the fluid 
description is justified, these can be given in the form 

r = r{D,v), u = u{D,v). 

In Section 3.1 we discuss, for illustration, the simplest possibility, adiabatic perturbations of an ideal 
fluid, where we just set 11 = F = 

There are however situations where the description of matter as a fluid is not sufficient. One 
then has to resort to the matter equations of more fundamental quantities, e.g. scalar fields and/or 
gauge fields. 



2.3 Collisonless Matter 

In this section we discuss another approximate description of matter which can be used for weakly 
interacting particles like photons in a recombined universe or massive neutrinos which might consti- 
tute the dark matter. Here, the basic quantity is the one-particle distribution function / which lives 
on the mass bundle, 

Pm = {{p,x) G TM\ g{x){p,p) = -m^} . 

When collisions can be neglected, the matter equation is the one particle Liouville equation (for a 
thorough treatment of the kinetic approach in general relativity see Stewart [1971] and references 
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therein). Choosing coordinates {x'^,p^) on (where p = p^di + p'^{x'^,p^)do), Liouville's equation 
reads 

{p^d, - r^^p-/±.)f = . 

In an unperturbed Priedmann universe, this equation is equivalent to the requirement that the 
distribution function, / is a function of the redshift corrected momentum, v := a^gijp^p^ alone. 

2.3.1 A gauge invariant variable for perturbations of the distribution function 

We want to split the distribution function / in a perturbed Friedmann universe into a background 
component and a perturbation. We cannot do this directly since the background distribution func- 
tion, / is not defined on but on the background mass bundle = {(p, x)| g(^{p,p) = —m?} . 
Therefore, to split /, we first have to choose an isomorphism i : P^n — ^ Pm '■ {x,p) i— > {x,p). Then we 
can define the perturbation according to 

f={f + F)ot. 

It is clear that the perturbation F in general depends on the isomorphism l which deviates in first 
order from identity. Choosing two basis (e^u) and (e^) which are tetrads with respect to g and g, we 
may set for p = p^d^j, = Tr'^e^ 

t(x,7r%) = {x,TT'^e^) . 

On the other hand, every isomorphism i is of this form (i.e. p{T^^e^) = ir^e^) for suitably chosen 
tetrads (e^) and (e^u). The tetrad e^^ can be defined by 

e,i = + R^e.^ , (2.65) 

where the symmetrical part of R is determined by the orthogonality condition: 

Rixv + Rvix = -a^h{ei^, e„) {g = g + a^h) . 

To determine how F transforms under gauge transformations, we consider a vector field X which 
defines a gauge transformation. The flow of X on A^, (f>^ , induces the flow T(f>^ on TM. The natural 
lift, TX, of X to TM is defined by T^f = (p^^ . A short calculation shows that for X = X^'d/^ 



TX = {X^d„X^^p^-^^) 

for a coordinate basis {xi^jp = p'^d^). For the full distribution function this leads to the transforma- 
tion law 

[/ : P„ ^ R] [{T^)J : r</)(P^) ^ R] . (2.66) 
In linearized form this yields 
/ ^ / + Lrxf . 

We now want to use our split of / to obtain a transformation law for F. The first problem to 

note here is that / is not defined on all of TA4 but only on PfYi: and since TX is in general not 
tangent to Pm, LtxJ is not well defined. But of course it is possible to extend the definition of / to 
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an open sub-set of TAi containing P^- We thus do not have to bother about this technical point. 
We just keep in mind that the gauge transformation properties of F should not depend on such an 
extension. 

More important is that F also depends on the choice of the isomorphism l, the transformation 
properties of which have also to be taken into account. We now choose two admissible splittings of 
/ given by 

f = flO Fi = f2 L2+ F2 . 

Then, there exists a gauge transformation given by a vector field X, such that /2 = /i — Lxxfi, and 
therefore 

fio ii+ Fi = fio L2+ F2- Lxxfi ■ 
The change of F under a gauge transformation is thus given by 

F2 - Fi = f o ii - f o L2 + LtxI = i{TX)|| Fi , 

where we have dropped the index 1 on /. 

If Li^2 are specified by the tetrads e(i^2)fi which are related to the background tetrad according 

to eq. ( p. 65 ) with matrices R^/}''^^'^ , we obtain 

/o .1 - fo C2 = ^XiR^^" - R^^n ^ -{TX)J , (2.67) 

d 

where we have introduced (TX)x = (i^jf^'' - R'^}'^")^^- 



For p = pf^d^ = n'^e^ we find 

(rx)ii =TX- (rx)x = x^'df, + x^.^p""-^ - (r^^> - R^^>)^i'— 



V 



dpi' ^ ^ ^ ' di^ 

It is easy to see that (rX)|j (^(p,p)) = which shows that (rX)|| is parallel to Pm-, i-e., (TX)||/ 
does not depend on the extension of / which is necessary to make {TX)f and (TX)±f well defined. 

We now explicitly calculate {TX)^^ for X = Tdt + L'^di and a tetrad (e^) which is adapted to the 
splitting of spacetime into {t = const.} hyper-surfaces , S^: 



eo = n , eo 



n = a ^dt , g{ei,n) = 



The triad (cj) is an orthogonal basis on the slices S^. We furthermore use the fact that there exist 
coordinates such that the metric of a space of constant curvature k is given by 

k 

7ij = (1 + -r'^y'^S.j = X^6ij . 

In this coordinates we can choose = {aX)~^di. According to our definition of the perturbed lapse 
function, shift vector and three metric, we then have 

eo = a-\dt-(3) = {l-A)eo-B'Xei 
Si = (1 - HL)ei - H- ^ej , 



31 



where the indices of the perturbation variables are, as usual raised and lowered with the metric 
7 and Hf is traceless, but may contain scalar, vector and tensor contributions. Using the gauge 
transformation properties of these variables and the definition of {TX)± we obtain 

(TXU = [(W + r]vr°A_A(rN-LV^ 

+[(d/a) - TS/ + 1/2(L,I^- + L\.)y^ . 
Inserting the expression for TX we find 



We now use that / is only a function of 



Denoting the direction cosines of the momentum by we obtain (e* = p'/ ^Y^i PiPi = ttVa/Z^i TTiTTj), 

op^ dv 

df idf 
-pr-^ = ae — 
OTT^ av 

Op^ OTT^ dv 

dp^ SttO ■ 
Furthermore, 

Thus, the terms containing U in {TX)\\f cancel. Introducing q = X^/v"^ + mP'o^ we finally obtain 

(rX)||/=Kd/a)r + ge^r,,]^. 

This leads to the transformation law 

F^F^ [v{d/a) + qe'di]T . (2.68) 
dv 
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We first note that F is invariant under vector type gauge transformations. Comparing ( |2.68| ) with 
the transformation properties of IZ and a given in ( ^.23| ) we find the following gauge invariant 
combination: 

F=F-\vn + qle'dia] ^ . (2.69) 

dv 

2.3.2 The perturbation of Liouville's equation 

Choosing an arbitrary basis of vector fields (e^) and corresponding momentum coordinates, p = vr^e^ 
on TA^, Liouville's equation is given by 

LxJ = TT'-e.if) - u^\ip)^^§^ = . (2.70) 
If we select, as above, a tetrad adapted to the slicing of Ai into slices of constant time, 

eo = n = a'^{dt - /3) , G TSt , 
we find (see Appendix A5) 

a 

_[^' (p _ + (^0)2^ _ ^-1(^1. _ cV)vrOvr^] A , 

-'a a J J (j^t 

where p = 7r*ej is the component of p tangent to the slices and c*^ is defined by 

The (i?*) are the basis of one forms dual to the vector fields (e*) on S^. More details are found in 
Appendix A. We now rewrite Xgf in terms of the variables 

t, X, V = a\/ TTjTT*, e* = an'^/v = and q = aXn^. For this we use the following easily established 
identities: 

(3 = B'di = aXB'ei 

a = a{l + A), a'^ = a'^{l-A) 

Since the background distribution function only depends on v, i.e., = oe*^, only the symmet- 
rical part of Pj'' contributes to the Liouville equation in first order and the term a;*j(/3) gives no 
contribution in this approximation. With the help of the splitting / = {f + F)ol and the background 
Liouville equation for / we finally obtain 

qdtF + ve'diF - vk\/2{x' - x^e^e')^ = [qh'A,i -qve'ejiB'^^ - if/) + v^Hl]^ . 
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For the gauge invariant tensor and vector contributions to F this yields, setting F^^) = jr(^) and 
qdtJ^iT) ^^^iQ^j:{T) j^kd^ = qve'eM^- (2.71) 

qd,T^y) +ve^d,TiV) -yT],eU'^^ = qve^e,a\P'^ . (2.72) 

For scalar perturbations the situation is somewhat more complicated. Since the scalar contribu- 
tion F^^^ to F is not gauge invariant we want to express the Liouville equation in terms of the gauge 
invariant combination 

dv 

After carefully calculating dt[vTZ + qle^a^i ], di[vlZ + qle^a^i ] and -^[vTZ + gZeV,j ], we finally obtain 
the Liouville equation for scalar perturbations in a Friedmann universe IPurrer, 199"0| ]: 

qd^jriS) ^ ^.ig.jriS) _ ^^^e^^^ = {q^d,^ - v^di^V^ ■ (2-73) 
2.3.3 Momentum Integrals 

To connect this equation of motion to Einstein's field equations, we calculate the energy momentum 
tensor from /, which is given by 

T''" = I P^p'fKp) , (2.74) 

where /i is an invariant measure on Pm{x) (for a general definition see Stewart [1971]). With respect 
to a tetrad p = ir'^e^, fj, looks like in special relativity 

= r-nl = -^dvdn , 2.75 



where we have used the definition of v and q, and is the usual surface element on the 2-sphere 
integrating over the momentum directions e. Let us, as an example, consider Tg : 



= p{l + 5) = - f ^-^{fiv) + F)dvdn 
J q 



Expressing p^ and po in terms of v and q and separating into a background and first order contribution 
yields 

47r /■ 2 - 1 /■ 2 

P = / V qfdv , p5 = / V qFdvdn . 

Using this we obtain 

I v'qJ^^^Uvdn = 6-^1 v^q^dvdn . 
pXa^ J pXa^ J dv 

After a partial integration, inserting ^ = )?v/q and 
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we end up with 

^ ( v^qT^^Uvdn = 5 + 3(1 + w)n = D„ . 



pXc 



(2.76) 



For the last equality sign, the definition of the gauge invariant density perturbation variable Dg is 
inserted. 

In a similar way all the other momentum integrals are obtained: 

-1 



lAV 



/2(n|,, - i/37,,An) 



a^p J q 



{l/3)5,j)J^^^Uvdn 



a^p J 3q A 



*p J 3q 

For the vector and tensor perturbations one finds 

-1 



(//2)(n 



(V) 



A '•-4 
^pj q 



-4- ! —{(-'e^ - {l/2,)5iAJ^^^Uvdn 
a^p J q 



A 

a^p 



(2.77) 
(2.78) 
(2.79) 

(2.80) 
(2.81) 
(2.82) 



These matter variables inserted in Einstein's equations ( 2.46| ) ( p. 48 ), ( 2.53| ) and (2.56) yield the 



geometrical perturbations <I>, cr* and H^J^ which enter in (|2.73| , 2.72 , 2.71 ). In Section 5, we discuss 



how this closed system is altered in the presence of seeds. 



2.3.4 The ultrarelativistic limit 

Here we briefly investigate the special case of extremely relativistic particles for which we can set 
m = 0. Since curvature only may play a role in the late, matter dominated stages of the universe, we 
neglect it here, A; = 0, A = 1, so that q = v. (The generalization to /c 7^ is straight forward.) In the 
extremely relativistic case all the integrals above contain the energy integral • / v^J-dvdQ,. Therefore, 
it makes sense to introduce the perturbation of the energy integrated distribution function, the 
brightness: 



pa* JO 



J'V^dv 



(2.83) 



which is a function of the momentum directions e* only. Defining 

47r r"" ^ o , 

Fv^dv , 



pa'^ Jo 



one finds, using (|]6^) and the gauge mvariance of F(^) and F(^) 
^{S) ^ i(^) + 47^ + 4/e*aiC7 



,{T) 



(2.84) 
(2.85) 
(2.86) 
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In the case where Ai describes thermal radiation, we may interpret the perturbation in the distri- 
bution function as a perturbation of the temperature: 

Inserting this form of / one obtains F = —v^ ■ The integral ( 2.83| ) then yields 

4 T 
4 T 

Therefore, (1/4) can be interpreted as a gauge invariant variable for the temperature perturba- 
tion]] 

In terms of Ai the perturbation equations (2.73, ^?72| , |2.71| ) become (for k = 0) 

AS) 



M^' + e'd,M'^^^ = 4e'di{<!>-'<i>) (2 



A^(^) + e'diM^''^ = -^e'ejalP' (2.89) 



M^^'^ + e'diM^^'^ = -Ae'e^Hij . (2.90) 

The evolution of the distribution of massless particles only depends on the Weyl part of the curvature. 
This is geometrically very reasonable since null geodesies are conformally invariant. 

By similar calculations like in the preceding paragraph, one finds the perturbations of the energy 
momentum tensor for extremely relativistic particles 

D = ^ f MdVL (2.91) 
47r J 

lAV = — I e'diMdn (2.92) 
167r J 

l^{n\,i - l/3e,,An) = ^ I {e'e^ - l6ii)Mdn (2.93) 

Att J 3 

r = (2.94) 

= — I e'M'^^Un (2.95) 
47r J 

(z/2)(n;;;)-n(;;)) = A|(ev-(i/3)5.,)A^(^)dn (2.96) 

njp = ^ ! {e'e^ - {l/3)6ij)M'^^^dvdn . (2.97) 



^Note that even though T, M and ST are scalar functions they do in general contain vector and tensor perturbations, 
since they depend not only on position but also on momentum or momentum direction. They may contain terms of 
the form o'e^ or r^-'eiej where a is a divergence free vector field and r is a traceless, divergencefree tensor field. These 
are the type of contributions which we indicate with and ^"^K 
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We use Liouville's equation for massless particles for the numerical calculation of perturbations of 
the cosmic microwave background in Chapter 5. In Section 3.2 we use ( 2.88 ) to derive the Boltzmann 
equation for photons in an electron proton plasma. 



2.4 The Propagation of Photons in a Perturbed Friedmann Uni- 
verse 

On their way from the last scattering surface into our antennas, the microwave photons travel through 
a perturbed Friedmann geometry. Thus, even if the photon temperature was completely uniform 
at the last scattering surface, we would receive it slightly perturbed [^achs and Wolfe, 1967 1. In 



addition, a photon traveling through a perturbed universe is in general deflected. In this section 
we calculate both these effects in first order perturbation theory. I present the calculation rather 
explicitly, since I haven't found a complete gauge invariant treatment of this problem anywhere in 
the literature. For sake of simplicity we restrict ourselves to A; = 0. 

As already mentioned, two metrics which are conformally equivalent, 

ds^ = a^ds^ , 

have the same lightlike geodesies, only the corresponding affine parameters are different. We may 
thus discuss the propagation of light in a perturbed Minkowski geometry. This simplifies things 
greatly. We denote the affine parameters by A and A respectively and the tangent vectors to the 
geodesic by 

dx dx 9 2 

n = —— and n = , n = n = , 
d\ dX 

with unperturbed values n*^ = 1 and = 1. If the tangent vector of the perturbed geodesic is 
(l,n) + 6n, the geodesic equation for the metric 

ds^ = {rjau + hay)dx^dx'^ 

yields to first order 

Jn'^lf = -r,'^-[/i,o + K,n']{ + ^ V,-^"^'*^^ > (2-98) 

where the integral is along the unperturbed photon trajectory and the unperturbed values for can 
be inserted. Starting from this general relation, let us first discuss the photon redshift. The ratio of 
the energy of a photon measured by some observer at to the energy emitted at ti is given by 

E,IE. = = (T,IT.)'f^ , (2^99) 

[n-u)i {n-u)i 

where uj and Uj are the four velocities of the observer and the emitter respectively and the factor 
Tf/Ti is the usual redshift which relates n and h. We write Tf /Ti and not aj/ai here, since also 
this redshift is slightly perturbed in general, and we want a to denote the unperturbed background 
expansion factor. 

Since this is a physical, intrinsically defined quantity it is independent of coordinates. It must thus 
be possible to write it in terms of gauge invariant variables. We now calculate the gauge invariant 
expression for Ef/Ei. The observer and emitter are comoving with the cosmic fluid. We have 

n = (1 - A)dt - Iv'di . 
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Furthermore, since the photon density p^^'^ oc may itself be perturbed 

Tf/T, = (a./a;)(l + ^ - ^) = a./a^l + (l/4)5(^)|f) , 
If li 

where 5^^^ is the intrinsic density perturbation in the radiation. This term was neglected in the orig- 
inal analysis of Sachs and Wolfe, but since it is gauge dependent, doing so violates gauge invariance. 
We therefore have to include 5^"^^ to obtain a gauge invariant expression. Inserting all this and ( 2.98 ) 
into (12^991 ) yields 

Ef/Ei = {ai/af)[l + n^v,j \{ + A\{ + (l/4)(5M|f - 1/2^^ ^nWA] . (2.100) 

With the help of equation ( 2.11| ) for the definition of h^y one finds for scalar perturbations after 
several integrations by part 

{Ej/E,)^^^ = {a,/aj){l + \{l/A)D^:^+lvlp\^ + nl- j\^-m>^} (2-101) 
= {ai/af){l + [{l/A)D^^'^ +ivlp^n^ + ^-^{- j^^{^-i>)d\} . (2.102) 

Here D^s \ Dg"^ denote the gauge invariant density perturbation in the radiation field and is the 
peculiar velocity of the matter component (the emitter and observer of radiation). From the second 
of these equations one sees explicitly that the geometrical part of the perturbation of the photon 
redshift depends on the Weyl curvature only (specialize eq. ( 2.26| ) to purely scalar perturbations), 
i.e., is conformally invariant. 

For a discussion of the Sachs- Wolfe effect alone we neglect the intrinsic density perturbation of the 

(r) 

radiation, i.e., we set Dg = 0, which now is a gauge invariant statement (but a bad approximation 
in many circumstances like, e.g. for adiabatic CDM perturbations). T/^™-) is a Doppler term due to 
the relative motion of emitter and receiver. The ^' — $ - term accounts for the redshift due to the 
difference of the gravitational field at the place of the emitter and receiver and the integral is a path 
dependent contribution to the redshift. 

For vector perturbations 6^'^^ and A vanish and eq. ( p.lO0| ) leads to 

{Ef/E,)(''^ = {a,/af)[l - v\-^\^i + a,n^d\\ . (2.103) 



Again we obtain a Doppler term and a gravitational contribution. For tensor perturbations, i.e. 
gravitational waves, only the gravitational part remains: 

{Ef/Ei)^^^ = {ai/af)[l - Hi^n^n^dX] . (2.104) 

Equations ( |2.101| ), ( |2.103| ) and ( p.l04| ) are the manifestly gauge invariant results for the Sachs- Wolfe 
effect for scalar vector and tensor perturbations. 

In addition to redshift, photons in a perturbed Friedmann universe also experiences deflection. 
We now calculate this effect. 

The direction of the light ray with respect to a comoving observer is given by the direction of the 
spacelike vector 



n 



(3) = n + {un)u 
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which hves on the sub-space of tangent space normal to u. Let us also define the vector field n^^^y 
which coincides with n(3) initially (i.e. at ti) and is Fermi transported along u, i.e. 



(nLv„ti)n 



(2.105) 



"(3) ~ ^'"(3) 

Note, that we have to require Fermi transport and not parallel transport since u is in general not 
a geodesic and therefore (nn^'g^ = is not conserved under parallel transport! (For an explanation 

of Fermi transport, see e.g. Straumann [1985]). n'^^ = (0, n) + Sn^j^^y where S-n)^^^ is determined by 
the Fermi transport equation ( 2.105 ). Since the observer Fermi transports her frame of reference 
with respect to which angles are measured, she would consider the light ray as not being deflected if 



n 



(3) 



(tf) is parallel to n(3)(tj). The difference between the direction of these two vectors is thus the 



light deflection: 



ipe 



n 



(3) 



K3) -^(3)) II 

II ^11 



(tf) 



(2.106) 



Here e is a spacelike unit vector normal to u and normal to n^^^ which determines the direction of the 
deflection and ip is the deflection angle. (Note that ( |2.106 ) is the general formula for light deflection 
in an arbitrary gravitational field. Up to this point we did not make any assumptions about the 
strength of the field.) For a spherically symmetric problem, as we shall encounter when discussing 
the collapsing texture, e is uniquely determined by the above conditions since the path of a light 
ray is confined to the plane normal to the angular momentum. In the general case, when angular 
momentum is not conserved e still has one degree of freedom. We now calculate ipe perturbatively. 
Let us recall and define the perturbed quantities: 



n 



u 



{l,n)+6n with = 1 



n 



(3) 



(1,0) +5u -- 
(0, n) + 5n(3') and 



1 + ^hoo,v) 



n 



(3) 



(0,n) +5n 



(3) 



The perturbation 5n is given in (2.98|). Furthermore, we obtain 



5n 



(3) 



with 



eu + 5n — 5u 



1 



(2.107) 



e = [n'v' - Sn" + -/iqo + n'hio] 



The Fermi transport equation leads to 

n^hio + Vi) 



^(nll/ 



6{n 



1 



(3)^ 



-[hijW 



+ 



/ 



dt{hjo,i - hioj)n 



(2.108) 
(2.109) 
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So that ('T'(3))^ = 1 and 

Lpe = (5n(3) — (nj-'g^ • n(3))n|g^ = eu + 5n — 5u — ^JT-j-g^ — [hjiU^n^ + n*((5(n|'gP* + 6n^ — v^)]n . 

Inserting ( p.l07l - ^.109| ) into (|2.106D we find 

cpe^ = (2.110) 
ipe' = 6' -{d- n)n' , with (2.111) 

dj = [6nj - Vj + ^hjkn'']\{ + ^ J,\hoj^k - hok,j)n''dt . (2.112) 

This quantitiy is observable and thus gauge invariant. For scalar perturbations one finds (after 
integrations by parts) 

i^,)^'^ = yj\{ + f{^-^),jdX. (2.113) 

For spherically symmetric perturbations, where e is uniquely defined, we can write this result in the 
form 

if = Vie'\{ + J^\<^> - ^),ie'dX . (2.114) 

The first term here denotes the special relativistic spherical aberration. The second term represents 
the gravitational light deflection. Here again one sees that gravitational light deflection, which of 
course is conformally invariant, only depends on the Weyl part of the curvature. As an easy test we 
insert the Schwarzschild weak field approximation: vj/ = — <i) = — ^i^. The unperturbed geodesic is 
given by X = (A,nA + ed), where d denotes the impact parameter of the photon. Inserting this into 
( |2.114 ) yields Einstein's well known result 

4GM 



For vector perturbations we obtain from ( 2.112 ) 

(5,)(^) = nj\{ - l[j\^j,k - (Tk,,Vdt + J^^ ak,jn^dX] . (2.115) 

This result can be expressed in three dimensional notation as follows: 
ipe = S — {S • n)n = 

-{U An) A n\{ + - f\'V A a) A ndt - /^(V(cr -n) An) A ndX . 

The first term is again a special relativistic "frame dragging" effect. The second term is the change 
of frame due to the gravitational field along the path of the observer and the third term gives the 
gravitational light deflection. This formula could be used to obtain in first order the light deflection 
in the vicinity of a rotating neutron star or a Kerr black hole. The special relativistic Thomas 
precession is not recovered with this formula since it is of order u^. 
For tensor perturbations we find 

(5,)(^) = -//jfcn^lf + f Hik,,rln^dX , 

Ji 
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or, after an integration by parts, 



(^e,)(^) = -Hjkn''\{ + / (Hikj + Hkinj)n^n''dX . (2.116) 

Ji 

Only the gravitational effects remain. The first contribution comes from the difference of the metric 
before and after the passage of the gravitational wave. Usually this term is negligible. The second 
term accumulates along the path of the photon. 



2.5 Gauge Invariant Perturbation Theory in the Presence of Seeds 

In this section we add an inhomogeneous term to the perturbation equations. Perturbations can 
then be generated even starting from an initially unperturbed spacetime. Seeds produce this in- 
homogeneous term in a natural way. By seeds we mean density perturbations originating from an 
inhomogeneously distributed form of energy whose mean density is much smaller than the density 
of the Friedmann background. We assume that, once they are produced, these seeds do not interact 
with the rest of the matter other than gravitationally. 



2.5.1 The energy momentum tensor 

Since the energy momentum tensor of the seeds, T/"j', has no homogeneous background contribution, 
it is gauge invariant by itself according to 

TI^^ can be calculated by solving the matter equations for the seeds in the Friedmann background 
geometry. (Since T^^ has no background component it satisfies the unperturbed matter and "con- 
servation" equations.) We again decompose T^X into scalar, vector and tensor contributions. They 
decouple within linear perturbation theory and it is thus possible to write the equations for each 
of these contributions separately. However, this decomposition is acausal. It requires T^jj[t,x) at a 
given time t to be known for all positions x and not only within a causally connected region. We 
just ignore this problem now and, nevertheless, work with this decomposition. We parametrize the 
scalar {S) vector iV) and tensor (T) contributions to T/^t' in the form 



= aV^^) = (MV/')/p (2.117) 

I^f) = aHv\f = {M^/l)f,\, (2.118) 

T,f^ = a^[(p(^)-(/V3)An(^))7., + /'ng] 

= m2[(/p//2_(i/3)A/,)7,,+/,|,j] (2.119) 

t/o''^) = a\f^ = {M^ /f)wf^ (2.120) 

Tf^ = ^[n£)+nW] = (MV20(^«![;.^+^?) (2.121) 

tP = a^ny = (MV/')r^^ . (2.122) 
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As before, / is introduced merely to keep the functions /., the vector fields w and the tensor field 
r"^ dimensionless. It may be chosen as a typical size of the seeds. M denotes a typical mass of the 
seeds. (It is of course possible to choose I = M~^.) 

(s) 

If we are given the full energy momentum tensor T^^ which may contain scalar, vector and tensor 
contributions, the scalar parts and /,r are in general determined by the non-local identities 

On the other hand A/„ and A(A + 3k) f.,^ are also determined in terms of fp and fp by the "conser- 
vation" equations: 

fp - lAf, + {a/a){fp + 3/p) = , (2.123) 
- Z/„ - 2(d/a)Z/„ + fp + (2/3)Z2(A + 3k)U = . (2.124) 
Once fv is known it is easy to get w'^ = P/M'^{TQi) — lfv,i- To obtain wj we use 

7j)l^_7j^)b = M!(A + 2fc)^f). 
Again w\ can also be obtained in terms of w\ with the help of the "conservation" euation: 

+ 2(^)^f ) - + 2fc)w;f ) = . (2.125) 

2.5.2 Perturbation equations 

The energy momentum tensor of the seeds is determined by the unperturbed equations of motion. 
The gravitational interaction with the perturbations of other components does not contribute to first 
order. We assume non-gravitational interactions with other components can be neglected. This is 
certainly a good approximation soon after the phase transition and thus can only affect the initial 
conditions. (Situations where non-gravitational interactions must not be neglected are discussed by 
Magueijo [1992].) The geometrical perturbations can then be separated into a part induced by the 
seeds and a part caused by the perturbations in the remaining matter components: 

^ ^ s ~l~ ^ rn ) ^ 

(s) , (m) 

TT, . _ rris) rr{m) 

Using Einstein's equations, we can calculate the geometry perturbations induced by the seeds: 

-{A + 3k)^s = e{fp/l^ + 3{a/a)U/l) (2.126) 
{a/a)^s-^s = efjl (2.127) 
A($, + *,) = -2eA/^ (2.128) 
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{a/a){4>s-[{l/a){^y]}+ 



{2a(d/o2)- + 3(d/a^)2}[^', - l/a( 



(A + 2k)a. 



2%r^ + ) 

a 



As) 



)■] = 2e{fp/l^ - {2/3)AU) 
- 2ewl^^/l^ 
-2tn}l^ll 



(2.129) 
(2.130) 
(2.131) 

(2.132) 



We assume that e = AtkGM'^ is much smaher than 1, so that linear perturbation analysis is justified. 
As before, the part of the geometrical perturbations induced by the matter are determined by equa- 
tions ( 2.46 ) to ( ^.48| ). But in the conservation equations and in any matter equations the fuh 
perturbations, Oi and B.ij have to be inserted. 

We now discuss the example of a single fluid with only scalar perturbations where 11 and V 
are given in terms of D and V . We assume that in addition to the seeds we have one perturbed 
matter component which we indicate by a subscript m. Other components which contribute to 
the background, but whose perturbations can be neglected, may also be present. The conservation 
equation ( 2.57D then reads 

Dra - 3wm{a/a)Dm = (A + 3A;)[(1 + Wm)lVm + 2{a/a)wml'^Ilm] - 3(1 + Wm)efv/l ■ (2.133) 
The last term describes the influence of the seeds. 



Solving this equation for (A + 3k)lVm. and inserting the result and its time derivative into ( 2.58 ) 
yields a second order equation for Dm- Using 

(A + 3k)^ = ATTGpma\Dm - 2wml\A + 3A;)n^) + e{fp/l^ + 3-fJl - 2(A + 3k)U) 



and the conservation equation (2.124) we flnd 



D-{A + 3k)clD + (1 + 3cl - 6w){d/a)D - {3w{d/a) - 9((i/a)2(c^ - w) + 

+{l + w)A7:Gpa^}D = 

(A + 3k)wr + 2{h/a)wl'^{A + 3k)U 

+{2{d/a)w - 6{d/af{ci - w) + {1 + w)8TTGa'^p + 2/3(A + 3k)w}l'^{A + 3k)U 

+{l + w)e{fp + 3fp)/l\ (2.134) 

where we have dropped the subscript m. 

This equation describes the behavior of density perturbations in the presence of seeds in an arbitrary 
Priedmann background. We have not used Friedmann's equations to express d/a in terms of w 
and d/a, or p in terms of (d/a)^ so that ( 2.134)) is valid also if there are unperturbed components 



which contribute to the expansion but not to the perturbation. Note that within this gauge invariant 
treatment the source term is, up to a factor {1+w), just the naively expected term 47rGa2(/j('^) + 3p*^*)) 
for all types of fluids. 
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We now simplify equation ( |2.134 ) in the case where 11 = F = (adiabatic perturbations and 



no anisotropic stresses) and k = 0: If one chooses a realistic density parameter 0.2 < l^o ^ 2, the 
curvature term can always be neglected at early times, e.g., for redshifts z > 5. It is of the order 
(max(/, Z/f )//fc)^ as compared to the other contributions. (Here / and Ik = k~^^'^ denote the typical 
size of the perturbation and the radius of curvature, respectively.) Under these assumptions, eq. 
( |2.134D becomes 



D - clAD + (1 + 3c^ - 6w){d/a)D- 
3[w{d/a)-3{a/af{cl-w) + {l+w){iTT/3)Gpa'^]D = S, (2.135) 

where S = {I + w)e{fp + 3fp)/P . 

We Fourier transform ( p.l35| ) (and denote the Fourier transform of D with the same letter): 

D + k'^c'iD + (1 + 3c^ - 6w){d/a)D 
-3[w{-d/a)-3{d/af{c^^-w) + {l + w){4TT/3)Gpa^]D = S. (2.136) 

S = (1 + w)e{ fp + 3fp)/P is the Fourier transform of S. 

From the homogeneous solutions Di and D2 of ( p.l36| ), we can find the perturbation induced by 
S with the Wronskian method: 

D = ciDi + C2D2 ] (2.137) 
ci = - [ {SD2/W)dt , C2= [ {SDi/W)dt , (2.138) 



where W = D1D2 — D1D2 is the Wronskian determinant of the homogeneous solution. 
This leads to the following general behavior: If the time dependence of Di, D2 and S can be 
approximated by power laws, D behaves like D oc 5 as long as 5 7^ 0. If Di and D2 are waves 
with approximately constant amplitude and frequency lo, D can be approximated by a wave with 
amplitude proportional to uj~^ J e^^^Sdt. Thus, only typical frequencies of the source finally survive. 
As a second example, we consider collisionless particles (again only scalar perturbations and 



A: = are considered). The source term on the r.h.s. of Liouville's equation ( 2.73 ) can be separated 
as above into a part due to the collisionless component and a part induced by the seeds. Equation 
(|2.73| ) then becomes 



{qdt + v''dk)r = '^[{q/v)v''du^^ - {v/q)v'dk^m] + S , (2.139) 



with 



S = ^[[q/v)v^dk^s - {v/q)v^dk^,] . (2.140) 
In the same way, one obtains for massless particles 

dtM + YdiM = Afdii^m -^m)+S , (2.141) 
with corresponding source term 

S = A^dii^s - *s) . 



44 



With the integrals for the fluid variables Dg, V, T and 11 given in Section 3 and Einstein's equations 
(|2.46| ) to ( 2.4^ ) for the geometrical perturbations ^'m, and induced by the collisionless component, 
this forms a closed system. 
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Chapter 3 



Some Applications of Cosmological 
Perturbation Theory 



3.1 Fluctuations of a Perfect Fluid 

A perfect fluid is free of anisotropic stresses, i.e. 11 = , Ilj = , Iljj = for scalar, vector 



and tensor perturbations, respectively. For scalar perturbations the dynamical equation ( 2.48 ) then 
relates the Bardeen potentials: 

^ = _$ . 

For vector perturbations we obtain from (2.54) 

a^cTj = const. , i.e. cjj = {a.^/a)^a* . 

The vector contribution to the shear of the equal time hyper-sufaces of perfect fluid perturbations 
thus decays like 1/a?. For tensor perturbations eq. ( 2.56| ) yields 



Hij + 2{h/a)Hij + {2k - A)Hij = , 

which describes a damped gravitational wave with damping scale d/a = t^^. 

The often used notion of isocurvature fluctuations is defined by ^' = <I> = ( for scalar pertur- 
bations), CTj = (for vector perturbations) and Hij = (for tensor perturbations) on super-horizon 
scales. Adiabatic fluctuations require F = 0. In a coupled baryon/photon universe this reduces to 
Dg^^ = {4/3)Dg~\ Note that fluctuations from topological defects are always isocurvature, since 
they emerge in a causal way from an initially homogeneous and isotropic Friedmann universe. 

Let us now solve the perturbation equations for adiabatic perturbations of a one component 
perfect fluid with w = =const. and negligible spatial curvature. This simplification is a good 
approximation during some periods of time (e.g. in the radiation dominated epoch). The very simple 
behavior of vector perturbations in this case is given by equation ( 2.64| ). For scalar perturbations 
equations ( |2.57| - [2l58 ) reduce to 



D-3w{h/a)D = {l + w)lAV (3.1) 



lV + (d/a)lV = — '—D + ^. (3.2) 

w + 1 



46 



Taking the Laplacian of ( p. 2] ) (using ( 2.4(^ ) and $ = — and inserting eq. (3J) we obtain a second 
order equation for D 



D + (a/a)(l - ?,w)D - [3?i;(a/a) + (1 + w)A'KGpa'^ - ciA]D = . 

For a spatially flat universe with —l/3<w = cl =const., the scale factor obeys a power law 

2 



(3.3) 



a = {i^VCt) 



with u 



and C = (87rG/3)pa2(^+i)/'^ = const. 



3w + l 

We now Fourier transform D, so that the Laplacian is replaced by a factor —k"^. Equation ( |3.3D can 
then be expressed as ordinary differential equation in the dimensionless variable rj = kt. The regime 
ry <C 1 describes perturbations with wavelength substantially larger than the size of the horizon and 
the regime rj ^ 1 describes perturbations with wavelength much smaller than the size of the horizon. 
Denoting the Fourier transform of D again with D we find 

{2 - iy){u - 1) - u{u + 1) 



,, 2z^-l , 



-D + wD = 



rj rj^ 
In terms of / = Dt]'^^'^ oc pa^D, this equation becomes 

rj rj"^ 

For w = ^ this is the well known Bessel differential equation whose general solution is 

/ = Aj^icsTj) + BnuicsTj) = Zy{csri) 
(see Abramowitz and Stegun [1970]). For our perturbation variables D,^ and V this yields 

D = rf-'' Z,{csrj) 



V 



V " Z^^i{csr])] , 



(3.4) 

(3.5) 

(3.6) 
(3.7) 

(3.8) 



where we have used Z^, = CsZy^\ — {v + \/ri)Zy [[Abramowitz and Stegun, 1970 1. This solution was 
originally obtained by Bardeen [1980]. The asymptotic behavior of Bessel functions yields 

= Ct]" + E?7-('^+i) , for Cs7] < 1 

and 

Zu = I cos(csr/ - a,,) + "I sm{csi] - a^) , for Csi] > 1 

with = 7r(i/ + l)/2. Therefore, the "growing" mode behaves like 



and 





= ^0 


D 


= D^rf 


V 


= ^or/ 




= ^oV- 


D 


= Dorj^ 


V 


= Vov^- 



for CgT] <C 1 



('^+^) COs(Cs?/ — tti. 

cos{csrj - Uy) 
cos{csrj - Uy^i] 



(3.9) 



> for CgT] » 1 



(3.10) 
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On scales below the sound horizon, CgT] ^ 1, density and velocity perturbations grow only if < 1, 
i.e. for w > 1/3 (or w < —1/2) and they decay for —1/2 < w < 1/3. Radiation {w = 1/3) represents 
the limiting case where the amplitude of density perturbations remains constant. Perturbations in 
the gravitational field always decay for > —1, i.e. for w > —2/3. On scales substantially larger than 
the sound horizon (cg?? <C 1) it might seem at first sight that density and velocity perturbations are 
growing. But one easily establishes that, e.g., the "growing" mode of the alternative gauge invariant 
density variable Dg = D — 3{l + w)'^ — 3{l + w){i'/r])V is constant. Therefore, in a coordinate system 
where Dg represents the density fluctuation (a slicing with TZ = 0), density perturbations do not 
grow. This shows that the behavior of density perturbations on these scales crucially depends on the 
coordinates and can not be inferred from the growth of the particular gauge invariant variable D. 
This leads us to define the perturbation amplitude A on super-horizon scales as the amplitude of the 
largest perturbation variable in a gauge where this quantity is a minimum. For scalar perturbations 
thus 

A= min (max{yl, i^T, 5, TTi, vry}) . 

{gauges} 

It is clear that this quantity is of the same order of magnitude as the largest gauge invariant variable. 
In our case therefore 

A^l^^l , 

i.e., super-horizon perturbations do not grow in amplitude, as one would also expect for causality 
reasons. 

Fortunately, this analysis, which shows that no perturbations with < tt; < 1/3 grow, does not 
hold for dust {w = = 0). In this case, equation (l3.3|) reduces to 

D+{2/r])D-{6/T]^)D = , 

with the general solution D = At^ + Bt^"^, yielding the growing mode solution 

D = DqT]'^ r] = kt 

y = Vbr/. 

Again on super-horizon scales, < 1 , the perturbation amplitude is given by ^'o and is thus constant, 
whereas on sub-horizon scales, 3> 1, the largest gauge invariant variable is = D^rj^, i.e., density 
perturbations do grow on sub-horizon scales. 

An additional important case is given by dust perturbations in a radiation dominated universe 
with pr ^ pm at times t <^ teq, where teq denotes the time when Pr = Pm, which happens at some 
time, because pr decays faster (oc a~^) than pm oc a~^. At t <^ teq the scale factor grows according 
to a oc t but A-nGprnO^ ~ 4iTGpra'^{t/teq) = (3/2)^. Equation 

D + il/t)D--^D = 0, 

which is approximately solved by 

D = Dologit/teq) , t<^teq. (3.11) 



then yields 
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This fact, that even dust perturbations cannot grow substantially in a radiation dominated universe, 



is called the Mezaros effect ||Mezaros, 1974 |. 



We can now draw the following conclusions: For adiabatic perturbations of a perfect fluid, in a 
universe where spatial curvature is negligible, the time evolution of D, ^ and V is given by 

^ f a , in a matter dominated background, , 

I log(t/t*) , for kt S> 1, in a radiation dominated background, 

^ = = const. . (3.13) 

For pure radiation perturbations {p = {l/3)p, F = 0, 11 = 0) one obtains 

j const., for A;t < 1 

I exp{ik{x - v/T73t)) , for /fct » 1 , ^ ' 

_ J const. , for kt <^ 1 Ci 

I ^exp(ifc(x-yT73t)) , forA;t>l. ^ ^ 

On super-horizon scales, kt <^ 1 the total perturbation amplitude is given by = |$|. Thus, even 
if one special variable like, e.g., D is growing the perturbation amplitude A as defined above remains 
constant. 

We note the important result: The only substantial growth of linear perturbation is that of 
sub-horizon sized, pressureless matter density perturbations in a matter dominated universe. Then 
D (X a and the gravitational potential is constant. 

3.2 The Perturbation of Boltzmann's Equation for Compton Scat- 
tering 

In this section, we restrict ourselves to scalar type perturbations, i.e., vector and tensor fields can be 
derived from scalar potentials. Furthermore, we set A; = 0. 

Following Section 2.3, the perturbed photon distribution function is denoted by / = {f + F) o l 
and lives on the one particle zero mass phase space, Vq = {{x,p) G TAi : g{x){p,p) = 0}. Choosing 
coordinates on Vq, Boltzmann's equation for / reads 



where C[f] is the collision integral depending on the cross section and angular dependence of the 
interactions considered. 

In terms of the gauge invariant perturbation variable A4 (see section 2.3.4) the coUisionless 
Boltzmann equation is ( |2.88| ) 

M + e'diM = Ae'dii^ - ^) > (3-16) 
with (1^ 



^yj- poo 

M = —r Tv^dv 
pa'^ Jo 

= t + 47^ + Ale'dia . 
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To this we have to add the colhsion term which is given by 

/ 3 J ^r^i / 3 J df+ df_ _ dt(e)+ di{e) 



pa^ J pa^ J dt dt dt dt ' 

where /+ and /_ denote the distribution of photons scattered into respectively out of the beam due 
to Compton scattering. 

In the matter (baryon/electron) rest frame, which we indicate by a prime, we know 



^{P,e) = ^Jf'ip',e'P{e,e')dn', 



where ng denotes the electron number density, or is the Thomson cross section, and lo is the nor- 
malized angular dependence of the Thomson cross section: 



/\2l ^ 3 , . , 1 



uj{e, e') = 3/4[l + (e • e') ] = 1 + ^^ij^j with = e^ej - -Sij . 
In the baryon rest frame which moves with four velocity u, the photon energy is given by 

p' = PljU^ ■ 

We denote by p the photon energy with respect to a tetrad adapted to the slicing of spacetime into 
{t = constant} hyper-surfaces: 

p = p^n^ , with n = a"^[(l - A)dt + ffdi] , see Chap. 2 . 

The lapse function and the shift vector of the slicing are given by a = a(l + A) and /3 = —B'^di . In 
first order, 

Po = ap{l + A)- apeiP^ , 

and in zeroth order, clearly. 

Pi = apsi . 

Furthermore, the baryon four velocity has the form = a~^{l — A) , = vPv'^ up to first order. 
This yields 

p' = p^,u^=p{l + ei{v'-fi')) . 

Using this identity and performing the integration over photon energies, we obtain 

The distribution of photons scattered out of the beam, has the well known form 
(see e.g. Lifshitz and Pitajewski [1983]) 

— = GTU^f {p ,e) , 
so that we finally obtain 
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where 5r = (l/47r) / i[e)dVt is the photon energy density perturbation. 

Using the definitions of the gauge-invariant variables M and V , we can write C in gauge-invariant 
form. 

C = arnelD^'^ -M + 4e^l^^V + ^eijU'^] , (3.17) 

with dJ''^ = (l/47r) / MdQ = Sr + 47^ and 
3 



M'^ = -J M{e')e',^dn'] . 



Since the term in square brackets of ( |3.17] ) is already first order we have C = = aC. So that 
the Boltzmann equation becomes 

M + e'diM = Ae'di{^ - ^) + aaxneiD^^'^ -M- AeHdiV + ^eijM'^] . (3.18) 

In the next two subsections we rewrite this equation for two special cases. 

Note that perturbations of the electron density, ne = ng + Sn-g do not contribute in first order 
for a homogeneous and isotropic background distribution of photons. The only first order term 
which accounts for the perturbations of baryons and therefore acts as a source term for the photon 
perturbations is the Doppler term 4e^diV which is due to the relative motion of baryons and photons. 



A comparison of the numerical integration of equation ( 3.18 ) with and without collision term (with 
a gravitational potential ^' — <1> originating from a collapsing texture) is shown in Fig. 7. There one 
sees that the collision term has two effects: 

1) Damping of the perturbations by several orders of magnitude. 

2) Broadening of the signal to about 7° (FWHM) which corresponds to the horizon scale at the 
time when Compton scattering 'freezes out', i.e. tx ~ t and the mean free path of photons becomes 
larger than the size of the horizon. In general, there will always be the question whether the damping 



term, D^^ — or the source term, Ae^diV, in equation (|3.18D wins. It seems that for standard 
CDM the source term is so strong, that no damping due to photon diffusion occurs [Efstathiou, 
private communication]. It is an interesting and partially unsolved problem, in which scenarios 
of structure formation, photon perturbations on small scales {9 < 7°) are effectively damped by 
reionization. It seems plausible to me, that models with topological defects, where perturbations are 
highly correlated, are affected more strongly than models with Gaussian perturbations. 

The collision term above also appears in the equation of motion of the baryons as a drag. The 
Thomson drag force is given by 

F, = / C[M]e,dn = ^^^(M, + mV) , (3.19) 

3 f 

with Mi = — e^Mdn . 



^ 47r 

This yields the following baryon equation of motion in an ionized plasma 

3pb 



WjV + {d/a)ldiV = di^ ^-^(Mj + AWiV) , (3.20) 



where we have added the drag force to eq. ( 2.58 ) with u; = = 0. 
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We now want to discuss equations (3.18, ^3o|) in the limit of very many collisions. Clearly the 
photon mean free path is given by tx = It = (acrT?i-e)~^- In lowest order tx/t and It /I these equations 
reduce to 



and 



= a-M,- = 3L» 



(r) 



(3.21) 



(3.22) 



where we made use of ( ^.21 ) and ( ^.28 ) below, for the last equal sign. Eq. ( ^.22 ) is equivalent to 
( |2.60| ) for radiation. Using also ( p.60[ ) for baryons, it; = 0, we obtain 

5 3 3 3 

This shows that entropy per baryon is conserved, F = 0. Inserting ( 3.21| ) in ( p.l8| ) we find up to first 
order in tx 



M = D^;^ -Ale'd.V + ^eijM'^ -tTlD^ -4k'diV + ^eijM'^ 



.g ^ ^l^J" I ^v. ^.^j^j-^ -^^ ^jy:^ . (3.23) 

Using ( |3.23D to calculate the drag force yields 

Fi = {pr/3)[miV - diD^;^ + 4dii<^ - ^)] • 

Inserting Fi in ( ^3oD , we obtain 

{pb + {A/3)pr)ld,V + pb{a/a)ldiV = (p,/3)9,Z?W + {pb + {4/3)pr)d.,^ - {4pr/3)di^ . 
This is equivalent to ( |2.61| ) iov p = p), + pr, p = Pr/3 and F = 11 = 0, if we use 

i^W = (4/3)Z?(™) and Dg = ^^'^^ + ^"^^^ . 

Pm ~l~ Pr 

In this limit therefore, baryons and photons behave like a single fluid with density p = Pr + Pm and 
pressure p = pr/3. 

From (|2^ and (|]6l|) we can derive a second order equation for Dg. To discuss the coupled 
matter radiation fluid we regard a plane wave D = D{t) exp(ifc • x). We then obtain 

b + clk'^D + (1 + 3c^ - 6w){h/a)b - 3[w{ii/a) - ((i/a)(3(c^ - w) - (1/2)(1 + w))]D = . 

For small wavelengths (which are required for the fluid approximation to be valid), l/tx ^ Csk ^ 
we may drop the term in square brackets. The ansatz D(t) = ^(t)exp(— z / kcgdi) then eliminates 
the terms of order c^fc^. For the terms of order Csk/t we obtain the equation 



2A/A + (1 - 3ci - 6w){d/a) + Cs/cs = . 



(3.24) 
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For the case = w =const. , this equation is solved by A oc {ktY with v = 2/{3w + 1), i.e., the 
short wave hmit In our situation we have 



w 



Pr 



3{pr + Pn 



d = , , = 4^1^ and 

^Pr+Pm) ^Pr + ^Pm 



Cs/cs = -3/2(a/a) 



Pv 



ipr + 3pm 

Using all this, one finds that 

p™ + (4/3K V^^ f p + pV^'' 



A 



Cs {Pr + PraYa'^ J \ CsP'^a^ 



solves ( 3.24 ) exactly, so that we finally obtain the approximate solution for the, tightly coupled 
matter radiation fluid 

D{t) oc (-^^) ' expi-ik fcsdt) . (3.25) 

Note that this short wave approximation is only valid in the limit t ^ l/(cs/c), thus the limit to 
the matter dominated regime, Cg ^ cannot be performed. In the limit to the radiation dominated 
regime, — > 1/3 and p oc we recover the acoustic waves with constant amplitude which we 
have already found in the last subsection. But also in this limit, we still need matter to ensure 

= l/{aaTne) <C t. In the oppotite case, t-r ^ t, radiation does not behave like an ideal fluid but 
it becomes collisionless and has to be treated with Liouville's equation ( p. 16 ). 

In the paragraph 3.2.3 we evolve M. up to order tj^ and obtain the damping of Ai by photon 
diffusion. 



3.2.1 Spherically symmetric Boltzmann equation 

In the spherically symmetric case M depends on the momentum direction e only via the variable 
H = {r • e)/r. Integrating out the independent angle if and choosing the coordinates in photon 
momentum space such that the third axis is parallel to r, one finds 



- / eijdip=\ , = (3,3) 

I (i,j) = (l,l)or(2,2) 



We thus obtain 



271" Q /"l 1 

'2 1/Q^/ ,''2 i/Q\A/(^'— /2 



dip[e,,M^^] = -ifi' - 1/3) / ifi" - l/3)Mdp' ^ -{p' - -)M2 
4 J-i 1 3 

For M{t,r,fj,) we further have 

1 - 

e'diM = pdrM + —d^M 
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so that Boltzmann's equation becomes 
1 - u2 

aarUelM - M- 4fidrV + (/x^ - 1/3) ^Mz] . (3.26) 

This is the equation which we have integrated numerically, (in somewhat different coordinates, see 
Chap. 4) for a $ — ^ from a spherically symmetric texture, to produce Fig. 7. 

3.2.2 Moment expansion 

As long as collisions are efficient, t > l/(ao"rne) = t^, it is reasonable to truncate M at second 
moments 

M = D^J'^ + CiM' + 5eijM'^ . (3.27) 
With this ansatz, the integration over directions of the zeroth, first and second moments of (3.18) 



yields 

+ ^diM' = , (3.28) 
+ djDl''^ + 2diM'^ = Adj{<^ - ^') - aaTne[M^ + AdjlV] , (3.29) 



M'^ + —{diM^ + djM') = -—aarrieM'^ . (3.30) 
15 10 



Taking the total divergence of ( |3?29| ) and ( ^^ , we obtain with the help of (|]2|) 



i)W - ^ADP -^a = A($ - ^) + aaTne[Dj[^ - ^AlV] (3.31) 

and 

Q - -ADJ^) = -—aarnea , (3.32) 
5 10 

where we have set a = didjM^^ . 



3.2.3 Damping by photon diffusion 

In this subsection we want to estimate the damping of CBR fluctuations in an ionized plasma using 
our gauge invariant approach, as it was done by Peebles [1980] within synchronous gauge. We again 
consider eqs. (|1|) and ( ^) , but since we are mainly interested in collisions which take place on 
time scales It ^ t, we neglect gravitational effects and the time dependence of the coefficients. We 
can then look for solutions of the form 

V (X Ai (X exp(i(fc • X — Lut)) . 
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In ( p.l8|) and ( 3.20D this yields (neglecting also the angular dependence of Compton scattering, i.e., 



the term eijM^^] 



M 



D 



(r) 



AikelV 



1 — itT{u} — k • e) 



and 



trkiolV = {pr/3pm)i4:iklV + M) 



(3.33) 



(3.34) 



with M = (3/47r) / eMdn. Integrating (|3.18| ) over angles, one obtains + {l/3)diM' = 0. With 
our ansatz therefore k • M = 3wDg^^ Using this after scalar multiplication of ( |3.34| ) with k, we find, 
setting R = 3pm/^Pr, 



IV 



LO — I 



A;2 ■ 



Inserting this result for V in ( 3.33 ) leads to 



1 , 

3 l-itT{oJ-kp) 



where we have set p = k • e/k. This is exactly the same result as in Peebles [1980], where this 
calculation is performed in synchronous gauge. Like in there (§92), one obtains in lowest order ujIt 
the dispersion relation 



R? + - ( R + 1) 

w = (Jo-i7 with Wo = A;/[3(l + i?)]^/2 ^nd 7 = (/c^t^/e) ^ 



{R+iy 



In the matter dominated regime, R^ 1, therefore 
7 PS k^tr/G . 



(3.35) 



(3.36) 



We now consider a texture which collapses in a matter dominated universe (a oc t^) at time 
tc = to/{l + Zc)^^"^. The total damping, exp(— /), which this perturbation experiences is given by the 
integral 



f « / 7it)dt . 

Jtc 



(3.37) 



The end time tend is the time, when the mean free path, tx equals tc, the size of the perturbation. 
We define Zdec as the redshift when photons and baryons decouple due to free streaming. This is 
about the time when the mean free path has grown up to the size of the horizon: txizdec) ~ t{zdec)- 
To obtain exponential damping (tc < t^^nd)-, we thus need Zc > z^ec- In this case damping is effective 
until 

l + Z,nd = {l+Zdecf'\l + Zcf'^ , 

and we obtain 



/ 



l + Zdec\^l^ 



1 + Zr 



1 + Zde 



15/8 



(3.38) 
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where we have set k'^ = (27r/tc)^- In terms of angles this yields 

f{e)^2{d/edfm/ey'^/''-i] , e<ed 

where 9d = + Zdec ~ 6°. 

If the plasma ionizes at a redshift Zj, z^ec < Zi < Zc, after the texture has already collapsed, 
damping is only effective after Zi, if Zi.nd < Zi. Instead of formula ( p. 38 ), we then obtain 



/ 



1 + Zd, 
l + Zc 



ec 



3/2 



1 + . X15/8 /i+ X15/8 



(3.39) 



This damping factor can again be converted in an angular damping factor in the above, obvious way. 
The CBR signal of textures is thus exponentially damped only if 

Z^ > Zdec ~ WOiOm/hoflBf/^ , 

and in this case only textures with 1 + Zdec < 1 + < (1 + -^j)^/ (1 + Zdec) are affected. So, if Zi < 50, 
there is little exponential damping, even if /150 = 2 and Qb = 0.1, but if Zi = 200, all textures which 
collapse at redshifts z^ > z^ are damped by a factor 10 or more. 

In this approximation, we have neglected the the amount of damping which still may occur after 
Zend and the induced fluctuations due to the source term oc IV in ( 3.1^ ). A more accurate numerical 
treatment where just equation ( 3.1^ ) is solved is shown in Fig. 7. 



Astonishingly, even perturbations on large scales up to the quadrupole may be damped by photon 
diffusion if the spectrum is steep enough (Peebles, private communication): We assume perturbations 
of a given size I are uncorrelated and have an average amplitude A. On a larger scale L = Nl, they 
statistically induce perturbations with typical amplitude AN~^^'^ = A{1 / L)^^"^ . (A cube of size 

contains A^^ cubes of size I. The statistical residual of perturbations with amplitude A is thus 

AAr-3/2.) 

This simple argument has two interesting conclusions: 

i) The effective spectrum \D{k)\'^ of Gaussian distributed fluctuations cannot decrease faster than 
|Z?(/c)p oc k"^ towards large scales. 

ii) In the limiting case, \D{k)\'^ oc k^, all the power in large scales is induced due to the statistical 
residuals of small scale perturbations. Therefore if small scale perturbations are damped so are 
large scale perturbations up to the size of the horizon. 

It is clear that this mechanism crucially depends on the assumption of uncorrelated perturbations. 
It is thus not effective for fluctuations induced by topological defects. It also does not work for a 
Harrison-Zel'dovich spectrum, \D(k)\'^ oc k. 



3.3 Perturbations of the Microwave Background 

The redshift of photons propagating in a perturbed Friedmann geometry is given by (for scalar 
perturbations) ( p.l02| ): 

^ = [iz)M - V,n= + - '^)dT . (3.40) 
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The first of tfiese terms is due to intrinsic fluctuations on the surface of last scattering, the second term 
is the usual, special relativistic Doppler shift, the third and the last terms are gravitational redshift 



contributions, the Sachs — Wolfe effect |Sachs and Wolfe, 1967|: The third term is the difference of 



the potential at the emitter and receiver and the last term is due to the time dependence of the 
gravitational potential along the path of the photon. 

Since perturbations of the cosmic microwave background (CMB) are probably the most reliable 
observational tool for investigating the initial perturbation spectrum, and since they are calculable 
within linear perturbation analysis, we present them here in some detail. There are seven different 
physical mechanisms which perturb the microwave background on different scales. The first four of 
them are given by equation ( |3.40| ): 

• Intrinsic inhomogeneities on the last scattering surface, 

AT 1 



r ) 



T 4 

• relative motions of emitter and observer. 



— = -V.n\l, 



• the difference of the gravitational potential at the position of emitter and observer 
(Sachs-Wolfe I), 

• and the time dependence of the gravitational field along the path of the photon 
(Sachs-Wolfe II), 

AT ■ 

— = (^ - $)dr . 

• In an intergalactic ionized plasma, fluctuations are damped by photon diffusion (Section 3.2.3). 
As long as the mean free path of photons is considerably smaller than the size of the horizon. 
It = l/{aTne) <C Ih, this damping increases exponentially with the damping rate 7 oc It/I'^ 
and is thus very effective for inhomogeneities with sizes smaller or of the order of the Thompson 
mean free path, I < 

,5T^ ,6T^ , ff 



(— )/ = (— )iexp(-y 7dt) 



Clearly, damping by photon diffusion can only be efficient as long as It < Ih- At 2; ~ 100, 
It ~ Ih (for ^2 = 1 and = 0.05) and angular scales which are larger than the horizon scale 
at z ~ 100, which corresponds to about 6°, are usually not affected (for an exception of this 
rule, see the note at the end of the last section). 
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On the other hand, the passage of photons through a cloud of ionized gas at a temperature 
Te different from the photon temperature T induces deviations in the black body spectrum, 
which can be cast in the so called Compton y parameter. This is the Sunyaev-Zel'dovich effect 
ISunyaev and Zel'dovich, 1973/80| ]. In the Rayleigh-Jeans regime the change of the spectrum 
corresponds to a temperature shift according to 



AT , ff.Te-T. 

-— w -2y , with y = / ( jUeaTdr . 

T Ji rUe 



• Since the Bardeen potentials of linear dust perturbations are constant, they do not give rise to a 
path dependent contribution to the photon redshift. But once a dust perturbation has become 
non-linear and is virialized so that its density pD ~ const. , the corresponding gravitational 
potential ^ grows like . Usually, the gravitational potential of such a dust cloud is still very 
small, so that the redshift of photons passing through it can be calculated in linear perturbation 
theory. Equation ( |3.40| ) then yields 

1 ■' 1 + z 

where Az denotes the redshift difference of the two ends of the structure and z is the average 
redshift of the perturbation. Typically not only ^ but also Az is very small. This non-linear 



contribution to AT jT is called the Rees-Sciama effect [ Rees and Sciama, 1968 1. 



From these seven effects, for a long time only the Doppler term, point two, was detected. It leads 
to the famous dipole anisotropy with amplitude {AT /T)dipoie ~ 2 x 10~^ which shows that we are 
moving with a relative velocity of about 600A;m/s with respect to the microwave background. Now, 
the COBE team confirmed that also the spectrum of the dipole anisotropy is the derivative of a 



blackbody spectrum with T = 2.TilK to an accuracy better than 1% [Mather et al., 199S] 
Recently, also the Sunyaev-Zel'dovich effect has been confirmed 



pcHardy et al., 1990| , [Birkinshaw et al., 1991]] , and with the COBE satellite [ [Wright et al., 1992 



pmoot et al., 1992 1, the sum of intrinsic fluctuations and the Sachs-Wolfe effect (points one, three 
and four) have been measured. With these recent observations, the cosmic microwave background 
has started to become a very successful tool for investigating cosmological perturbations on the linear 
level. 



3.4 Light Deflection 

We now want to present some applications of the effect of light deflection in a perturbed Friedmann 
universe. The general formulae are derived in Section 2.4. 

3.4.1 Monopoles 

As a first example, we discuss light deflection and lensing in the field of a global monopole, see also 
Barriola and Vilenkin [1989]. We discuss the simple, static hedgehog solution of a three component 
scalar fleld with 0^ = r/^, i.e., a non-linear sigma model on S'^: 

= ijx^ /r . 
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This is an infinite action and infinite energy solution and should thus not be taken seriously at large 
distances. In a cosmological context, when monopoles form via the Kibble mechanism during a 
symmetry breaking phase transition, the hedgehog solution may be approximately valid on distances 
small compared to the distance to the next monopole or antimonopole, which is about horizon 
distance. This is also the scale where the approximation of a static, i.e., non-expanding background, 
which we adopt here, breaks down. 

The energy momentum tensor of the hedgehog solution is readily calculated: 



too 



2 / 2 

7] jr 

Setting rf = we obtain, using the definitions ( 2.117 ) to ( 2.119 ) 



/p = l^r^ (3.41) 

A = (3.42) 

/p = -(l/3)/p (3.43) 

U = -(l/2)log(r/0 . (3.44) 

It is interesting to note that the quantity fp + 3/p which enters as source term in the evolution 
equation for density perturbations ( 2.134| ) vanishes, which shows that static global monopoles do not 



produce an attractive gravitational force, much like cosmic strings. 
Setting e = SvrGr?^ we find from (|2T26D and (|2T28D 



$ = -elog(r//) , ^ = . (3.45) 

The second of these equations again shows that the analog to the Newtonian potential vanishes for 
global monopoles, like it does for cosmic strings (exercise for the reader). The divergence of $ at 
large distances reflects the inflnite energy of our solution which needs a physical cutoff (at most at 
the distance to the next monopole). 

We consider a photon passing the monopole at t = 0, say with an impact parameter h. Its 
unperturbed trajectory is then given hy x = An + be. According to ( p.ll4| ), neglecting spherical 
aberration, we obtain 

ip = + - ^) ie^dX = -e TnT^dX = -e arctan(A/6) -evr = -Svr^Gr?^ . (3.46) 

Ji ' + 

This result was originally obtained (by completely different means) by Barriola and Vilenkin [1989]. 

We want to investigate the situation of gravitational lensing. First, we treat the special case, 
where source (S), monopole (M) and observer (O) are perfectly aligned at distances SM = s and 
MO = d from each other. If emitted at a small angle a = b/ s, a photon will reach the observer at 
an angle (3 = b/d, if the deflection angle \ip\ = a + (3 = b{s + d)/sd. This leads to an Einstein ring 
with opening angle 

f3 = b/d= \ip\s/{s + d)= M'^Gjfs/{s + d) . (3.47) 

If observer and source are slightly misaligned by an angle less than |(/?|, the ring is reduced to two 
points with the same angular separation. 

For monopoles produced at a typical GUT scale, rj lO^^GeF, the deflection angle is |(^| 10 
arcsec and thus observable. Since the density of global monopoles is about one per horizon citeBR, 
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roughly 10 monopoles present at a redshift z = A would be visible for us today. The probability that 
one of them is within less than 10 arcsec of a quasar with redshift z > 4, so that the lensing event 
discussed above could occur, is very small indeed. 



3.4.2 Light deflection due to gravitational waves 

We discuss formula ( |2.116| ) for light deflection due to a passing gravitational wave pulse for which 
the difference of the gravitational field before and after the passage of the wave is negligible: 

V^ej = J,\Hik,j +Hiknj)n^iT}d\ . (3.48) 

We consider a plane wave, 

Hki = ^{eki exp(i(fc • x — ujt)) with ekik^ = . 
For a photon with unperturbed trajectory x = Xq + An , t = A, we obtain 



^ieiik ■ a;o)e;^^/^"._fc2z:^^i_ [g(*(fc • n-^)Xf) _ ^{i{k ■ n^M]] for k ^ ton 
for fc = ton 



Setting n = {p/uj)k + qn± with = 1 and + = 1, we have eim'nJ'n'^ = q'^£±, where e± = 
eijn\n'^. Inserting this above, we find for the deflection angle an equation of the form 

/2„2 

V? = e^-^J=cos(a + w(p- l)t) . (3.49) 

Here e_\_ is determined by the amplitude of the gravitational wave and q , = 1 — (f' \s determined 
by the intersection angle of the photon with the gravitational wave as explained above. 

This effect for a gravitational wave from two coalescing black holes would be quite remarkable: 
Since for this (most prominent) event e± can be as large as ~ 0.1(-Rs/r), the rays of sources behind 
the black hole with impact parameters up to 6 < lO^i?^ would be deflected by a measurable amount: 

VP ^ 2"{W'Rs/b) . 

Setting the the source at distance dLS from the coalescing black holes and at distance ds from us, 
we observe a deviation angle 

a = (pdLs/ds ■ 

The best source candidates would thus be quasars for which d^s/ds is of order unity for all coalescing 
black holes with, say z < 1. In the vicinity of the black holes (r < WRs), linear perturbation theory is 
of course not applicable and also the 'monopole' component of the gravitational field is not negligible. 
But in the wide range W'^Rs > b > WRs (for radio sources) and W^Rs > b > WRs (for optical 
sources) our calculation is valid and leads to an effect that is in principle detectable. 

A thorough investigation of the possibility of detecting gravitational waves of coalescing black 
holes out to cosmological distances by this effect may be worth while. 
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Chapter 4 

Textures in Flat Space 



4.1 The (7— Model Approximation for Texture Dynamics 

Global texture occurs in every symmetry breaking phase transition where a global symmetry group 
G is broken to a sub-group H, such that 773(0 /H) / 1. 

In the cooler, broken symmetry phase the Higgs field remains most of the time in the vacuum 
manifold M = G/H, where the effective potential assumes its minimum value. It only leaves M if 
the gradients, i.e., the kinetic energy of the massless Goldstone modes, become comparable to the 
symmetry breaking scale ry. In order to discuss the dynamics of the massless modes at temperatures 
T <C Tc r/, it is thus sufficient to neglect the potential and, instead, fix the Higgs field in M with a 
Lagrange multiplier. 

For illustration, and since we believe that the results remains valid at least qualitatively also for 
other symmetry groups, we discuss the simplest version, an SU (2) symmetry which is completely 
broken by a valued scalar field <p. In this case, clearly M = SU{2) = S^. Here S" denotes the 
n-sphere and = means topological equivalence. 

The Higgs field (f) is described by the zero temperature Lagrangian 

L=^cP,^^'''-X{4>-4>-r]^? , (4.1) 

with corresponding field equation 

n(f> + AX(l){(f> ■ ^ - r}"^) = . (4.2) 

As argued above, on energy scales well below the symmetry breaking scale rj we may fix (fx^ = by 
the Lagrangian multiplier A((?!)- (j) — rj'^). Instead of the usual equations of motion, we then obtain the 
scale invariant equation 

n(f>-n-^{(l)-n(t))(f> = o , (4.3) 

i.e., (f) describes a harmonic map from spacetime to S^. If ^ is asymptotically parallel, (f>{r) -^r^oo 4>0j 
we can regard it (at fixed time) as a map from compactified three dimensional space (= S'^) into S'^. 
The degree (winding number) of this map is called the texture winding number. It is given by the 
integral of the closed three-form 



61 



{lo is nothing else than the pullback of the volume form on S'^). This integral over some region of space 
(e.g. a horizon volume) is of course well defined also if (p is not parallel at infinity and is often referred 
to as fractional texture winding number. Numerical simulations show that a texture starts collapsing 
as soon as the fractional winding inside the horizon exceeds about 0.5 [Leese and Prokopec, 1991 



Borill et al., 1991| ]. 

In our case even an exact spherically symmetric solution to equation (^^) is known: A spherically 
symmetric ansatz for (f) is 

(/>(r, t) = r]{r sin x, cos x) , (4.4) 

where r denotes the unit vector in direction of r and x is ^in angular variable depending on r and t 
only. In flat space the field equations for (j) yield 

- dh + dh + Idrx = ^ . (4.5) 
We now look for solutions which depend only on the self similarity variable y = {t — tc)/r, where 



tc is an arbitrary time constant. In terms of y, equation (4.5) becomes 
(y2_l)x" = sin2x , 

with the exact solutions x = 2arctan(iby) it rnr which were found by Turok and Spergel [1990]. To 
describe a collapsing texture which has winding number 1 for t < tc and winding number for t > tc, 
we patch together these solutions in the following way: 

, , _ J 2 arctan y + vr , — oo < y < 1 {a ei\ 

" I 2arctan(l/y) +7r , 1 < y < oo . 

It is straight forward to calculate the integral of the density to given above for this example and one 
of course obtains 



u{t) 



1, if t < tc 
0, if t > tc 



The kink of x at y = 1 reflects the singularity of the cr-model approach at the unwinding event 
t — tc = r = 0. There the gradient energy of the solution ( [4.6D diverges, i.e., becomes bigger than 
the symmetry breaking scale rj. Therefore, the Higgs field leaves the vacuum manifold, unwinds and 
the kinetic energy can be radiated away in massless Goldstone modes. To remove the singularity at 
r = 0, t = tc, we would have to evolve the innermost region of the texture (r < r/^^) with the true 
field equation (|4.2| ) during the collapse {tc — ri~^ < t < tc + 'r]~^). It is well known [Derrick, 1964 , 



Derrick's theorem], that a cr-model like ( [4.3| ) cannot have static finite energy solutions, but it is 
not proven that there is no static infinite energy solution. It is easy to see that (|4.6| ) represents 
an infinite energy solution. We shall see that the infinite energy of the solution will require some 
'renormalization' or, equivalently, a physical cutoff. 

Remark: It seems quite natural that this scale invariant problem admits scale invariant solutions. 
But this is by no means guarantied. It is well known, e.g., that Yang Mills theories in general do not 
admit non-trivial spherically symmetric solutions. But also other cr-type models do not admit them. 
As an example we have studied SO (3) broken to SO {2) by a vector field (the Heisenberg model). 
The vacuum manifold in this example is with 7r3(S^) = Z (the Hopf fibration). A spherically 
symmetric ansatz for a texture configuration in this model is the Hopf map: 

7i :R3^S2:r^n = i?(0(r),r)no , (4.7) 
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where R{0, e) denotes a rotation around e with angle 9, r is the unit vector in direction of r and rio 
is an arbitrary but fixed unit vector. For the mapping ( [4.7D to be well defined, we must require the 
boundary conditions 

6»(r = 0) = , 9{r = oo) = NItt . (4.8) 
is the Hopf invariant, or the ir^, winding number of Ti. The CT-model equations for n are 
□ n - (n • Dnjii = 0. (4.9) 
With the ansatz ( [4.7D above, this gives rise to the following two equations for 6: 

□0-^ = (4.10) 

-id,9)' + {dref-^^^-^ = 0, (4.11) 

which have no common non-trivial solutions. 

4.2 Gravitational Effects of Textures in Flat Space 

As long as e = IQirGrf' is much smaller than 1, the gravitational field of a texture is weak and we can 
calculate it in first order perturbation theory. We apply the gauge invariant formalism developed in 
Chapter 2. 

The energy momentum tensor, 

= 0,. -l/2g^,(t),x 4>'^ , (4.12) 
of solution ( [4. 61 ) is readily calculated with the result 
_ 2r/^ l + 3y' 



(1 + 
Arf y 



r2 (l+y2)2 

2rf 1 - y2 
(1 + y2 



(4.14) 



= ^IT^T^^'r (^-^^^ 



Clearly ( [4.6D represents a solution with infinite action and infinite energy, 
rR 



E{R) = 4tt / Toor^dr oc R , for R » t - tc 
Jo 



diverges linearly. 

Since this texture solution is spherically symmetric, it only gives rise to scalar perturbations. We 
now set = Arj'^ and I = tc, the 'radius' of the texture (for an expanding universe, tc is the horizon 
scale at texture collapse). We then find 

_ l + 3y^ 

JP - 2r2 (1 + 2/2)2 l^-^'^^ 
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fv = -^-^ (4.17) 
f 1 - y2 

= 0. (4.19) 

Prom (4.16) to ( [4.19| ) and the equations for the gauge invariant Bardeen potentials of scalar seed 
perturbations ( |2l26| , ^l2^ , [2l28| ), we obtain in flat space (a = , a = 1) 



e 1 + 37/2 
2r2(l + y2)2 



2 r(l+y2) 
A($, + ^,) = 0, 
with the solution 

= -\H{l + y\^/tl) (4.20) 
e 1 + 

M/, = jln(^). (4.21) 
^'s is only determined up to a function of time, which we have chosen to ensure ^ s — > 0, for t ±00. 

Of course, physical observables do not depend on this choice. 

It might seem unphysical that the potentials $s and ^ s do not vanish at infinity r — > 00, in 
contrary, they diverge. This reflects again the infinite energy of our solution. Noting this divergence, 
one might fear that linear perturbation theory breaks down at large distances from the texture, but 
we find that the relevant geometrical quantities, like e.g. the 3-dimensional Riemann scalar on the 
surfaces of constant time, 



5^R = 4fl-^A7^ ^ 4a-^A<D = + ^.^^ , (4.23) 

(r^ + [t — tc)^ 



do vanish at infinity. So that, far away from the collapsing texture and at early and late times, 
the solution does approach fiat space. The validity of linear perturbation theory for textures is also 
confirmed in Durrer et al. [1991], where we find an exact solution for a texture coupled to gravity 
and show that, for e < 0.1 say, it deviates only very little from the flat space solution used here. 

We now calculate the behavior of baryons and collisionless particles (dark matter or photons) 
and photons in this geometry. 
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4.2.1 Baryons around a collapsing texture 

Let us briefly discuss the behavior of cosmic dust (baryons) in the field of a texture. Equation ( p.l34| ) 
for a flat dust universe (c^ = u) = , d = 0) yields 

^ + Ai^GptlD = S , (4.24) 

with 

S = 2e- ^ 



'(l + y2)2 • 

The term AirGpt^D is the coupling of the perturbation to its own gravitational field. It leads to 
exponential growth of perturbations which is a feature of the non-expanding universe only. But our 
approximation, neglecting expansion, means that all times involved are much smaller than Hubble 
time. This coincides with AirGpt'^ <C 1. Within our approximation, it is thus consistent to neglect 
the self gravitating term in ( [4.241) . Direct integration then yields the solution 

D = e[yarctg(y) + ciy + C2] 

= e(t/r)[arctg(t/r) + 7r/2] + e , (4.25) 

where we have chosen the integration constants ci and C2 such that D converges to for large negative 
times, D(t = —00) = 0, and D converges to a constant for large radii, D(t, r = 00) = e. Since D 
is only a function of the self similarity variable t/r, we cannot consistently choose both boundary 
conditions to be 0. For late times, t/r ^ I, D grows linearly with time: 

D = eTT{t/r) , 

Near the time of collapse, \t/r\ <^ 1, D is of the order of e, D{t = 0, r) = e. A given time after 
texture collapse which is small compared to the Hubble time, D has the following profile: For large 
radii D k, e and roughly at r = bends into D ~ ei^t^jr and diverges for r — > 0. This divergence 
leads to early formation of non-linear structure on small scales. At time perturbations on scales 
of the order of r < r„j = eTrt^, have become non-linear. 

The total mass accumulated around a texture diverges like the mass of the texture itself (see [4]). 
But in the real, expanding universe one has to cut it of at roughly the Hubble radius at the time 
when the texture collapses. In- 

In this simple approximation, we end up with the following picture: Due to textures forming at a 
time t in the universe, objects of mass M ~ 2eM//(t), form at separations on the order of p~^l^lu{t\ 
Where Mu denotes the horizon mass at the time when the texture collapses and p is the probability 
that a four component vector field which is distributed in a completely uncorrelated manner over 
a 2-sphere winds around a 3-sphere (i.e. the probability of texture formation at the horizon). In 
Gooding et al. [1991], Aminneborg [1992] this probability has been found numerically to be about 
1/25. 



From matter conservation, (2.133) for w = we obtain 
AF = e/r[arctg(t/r) + 7r/2] - (e/2)- * 



r2 + 12 
and therefore 

Vi = -diV = -i!l[arctg(t/r) + tt/2]. (4.26) 
2 r 
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The total change in a particle's velocity as the texture collapses is thus independent of the particles 
distance from it and is given by 



/^Vj = Vj{oo) — Vj{ — oo) 



r 



(4.27) 



A result which was found by Turok and Spergel [1990] and Durrer [1990]. 

A numerical calculation for the distribution of texture in an expanding Friedmann universe, where 
the growth of density perturbations is given according to ( 2.135 ) with w = = 0, is presented in 
Gooding et al. [1991]. 



4.2.2 Collisionless particles in the gravitational field of a texture 



Let us now calculate the perturbations in the distribution function of collisionless particles induced 
by a collapsing texture. Like for dust, we neglect self gravity. We start with the gauge invariant 
perturbation equation for Liouville's equation in a Priedmann universe with k = (|2l39|) : 



qdtJ' + v^dkJ'=^[{qyv)v^dk^s-vv^dk'^s]=S 
av 



(4.28) 



where and $<j are the metric perturbations due to the texture. 

Making use of spherical symmetry and inserting the results ( ^.22| ) and ( [4.20 ) for 4' and $ yields 



qdtT + vndrJ" + ^ —df,J^ = S 



(4.29) 



with 



df 



dv r(l + y^) 



The solution of ( [4.29| ) is easily found with the help of the physical coordinates 
T = t—{q/v)fir , b = r\/l — fi'^ and t. 



(4.30) 

It is straight forward to see that r is just the time of closest encounter of the given particle with the 
texture, the impact time, and b is the impact parameter. In these coordinates eq. ( 4.29| ) reduces to 



qdtJ^{t, T,b) = S{t, T, b) and thus 

1 /■* 
Ht,r,b) = - 

q Jto 

A short calculation gives 



5(t,r, b)dt 



(4.31) 



Tqv 



^/v^T^~+Jq^Tv'^ 



= arctan 



h + v'^{t-T) 



q\Jv^T^ + (g^ + v'^yS^ 



(4.32) 
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We are interested in the total change of the distribution function due to the collapsing texture, i.e. 
in the limit ti — > — oo and t2 
this limit and the difference, 




oo. The logarithmic term in eq. ( 4.32 ), let us call it L, diverges in 



crucially depends on how we perform the limit. It can take any value from — oo to oo. This is due to 
the fact that we are dealing with an infinite energy solution, and we certainly have to 'renormalize' 
our results. If we would change the energy momentum tensor in a way that the texture would be 
'born' some time in the finite past, or if we would compensate it in a consistent way, as we do it in the 
expanding universe, this problem would disappear. A physically intuitive procedure is to introduce 
a cutoff at some time t2 = —ti ^> b, |r|. With such a cutoff the logarithmic term cancels, and we 
obtain 

e df Tqv _ ^ / q'^t + v'^ {t - t) 



J='{t,T,b,v) 



2 dv y/v^T^ + (g2 + t;2)52 



arctan 



We then can remove the cutoff and obtain the change of the distribution function long after the 
corresponding particles have passed the texture 



T{t,T, b,v) 




Tqv 



y/vh^^^{q^ + W)W 



(4.33) 



Introducing our old variables 

r2 = (t-rfv^q' + b^ 



we find 



it 



T 



qr 



J='{t,r,^,v) = — - — 




q{vt — qfir) 



. (4.34) 

■\/ {vt — qfir)'^ + r2(l — pi^){q^ + v'^) 

This result can be inserted in equations ( |2.76 ) to ( p.78| ) to obtain the induced perturbations of the 
energy momentum tensor. Here we just perform the non-relativistic and extremely relativistic limits. 



Non relativistic limit: 

In the non-relativistic case we have v <^ q. The gravitational field of the texture is strong only for 
r ~ \t\. We can therefore expand (|4.34 ) in the small quantity vt/{qr). The first order approximation 
yields 

T{t,r,fi,v) = {eTT/2)^^[^,+ {vt/qr){fi^ -l) + 0{vt/qrf] , (4.35) 
This leads to 

V = — [ J^v^fidvdfi = -— (4.36) 
p J 2 

27]- r 

D = — J^v'^qdvdp = eixt/r , (4.37) 
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the well known late time {t ^ r) results of the preceding subsection. 



Mit, r, /.) = 2evr^===^__ . (4.39) 



Extremely relativistic limit: 

In this case we have v = q and therefore 

Integrating J^v"^ over v and dividing by p/An yields the fractional perturbation of the brightness 
introduced in Chapter 2, 

t — fir 

Since M describes the fluctuations in the energy density radiated in a given direction /i, the 
temperature fluctuation is given by AT/T = (see also the more explicit discussion of this 

point below eqn ( |2.83| )). This coincides exactly with eq. ( [4.43 ) for photon redshift in the next 
paragraph. 

In a flat, eternal universe the signal from a texture collapsing at f = 0, as seen from an observer 
at time to and distance Tq would thus be 

^(0) = f , , . («0) 

A similar calculation using a specific gauge is presented in Durrer et al. [1992a]. Unlike in the 
expanding universe (see next chapter), there is no horizon present in this calculations. Photons that 
pass the texture long before or after collapse, \to\ ^ Vq are still influenced by it and yield even a 
maximum temperature shift, 
AT _ _^7r 

~Y ~ 2^ ■ 

In the expanding universe of finite age we expect, because of the finite size of the event horizon, 
AT/T to achieve a maximum for ~ to — tc and to vanish for to — tc <^ Tq. The comparison of the 
flat space result ( fl.40 ) and the effect of a compensated texture in the expanding universe is shown 
in Fig. 8. 



4.2.3 Redshift of photons in the texture metric 



Let us first determine the energy shift which a photon experiences by passing a texture. Without loss 
of generality, we set tc = in this paragraph. If we neglect the distinctive dipole term and intrinsic 
density perturbations, equation ( p. 101 ) leads to 



6£ 



(6 - <^)d\ + ^ 



(4.41) 



Denoting the impact parameter of the photon trajectory by b and the time when the photon passes 
the texture (the impact time) by r, we get r'^ = 6^ + (t — r)^ . Eq. (4.41) then yields after the same 
"renormalization procedure which led to (4.33) 



5^ 

'e' 



f 



2(r2 + 262)1/2 



[arctg( 



2t 



Vr2 + 262 



)] 



/ 



(4.42) 
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For tf, —ti » T, b, we obtain 



6£ 



(W2) 



(r2 + 262)1/2 



(4.43) 



This result was first found by different methods by Turok and Spergel [1990]. Photons which pass 
the texture before it cohapses, r < 0, are redshifted, and photons passing it after cohapse, r > 0, 
are blueshifted. This produces a very distinctive hot spot — cold spot signal in the microwave sky 
wherever a texture has collapsed. 

Of course our result is not strictly correct in the expanding universe, since we have neglected 
expansion in the calculation of and But the main contribution to the energy shift comes 
from times \t\ < \t\ + b. Therefore, our approximation is reasonable also for the expanding case, if 
|r| < Ih and b < Ih, where Ih denotes the horizon distance at the time of collapse, t = 0. On the 
other hand, by causality the texture cannot have a big effect on photon trajectories with |r| > Ih or 
b > Ih- a first approximation to the situation in the expanding universe is thus 



6E_ 

'e' 



f 



Vr2+262) 



for |r| < Ih and b < Ih 
for |r| > Ih or b > Ih ■ 



(4.44) 



4.2.4 Light Deflection 

To obtain the light deflection in the gravitational field of a spherically symmetric texture, we consider 
(as above) a photon passing the texture at impact time r at a distance b (impact parameter). The 
trajectory of this photon is then given by x(A) = (r + A, An + be). Making use of eq. ( |2.114| ), 
neglecting spherical abberation due to the relative motion of emitter and observer, we find 



<P 



Ji 



62 + 2A2 + 2Ar + r2'^^ 



CO 



= eTT , ^ . (4.45) 
V262 + r2 ^ ^ 

For T = 0, the deflection angle assumes the maximum value, v^max = e7r/\/2. For 6 = or r ^ oo. 



light deflection vanishes. This result was first obtained by different methods in [purrer et al., 1992a |. 



In contrast to global monopoles and strings, the deflection angle is not independent of the impact 
parameter b of the photon trajectory (except for photons which pass the texture exactly at collapse 
time, r = 0). This leads to a qualitative difference to the lensing caused by global monopoles and 
cosmic strings [ Barriola and Vilenkin, 1989 , [Vilenkin, 1984 ]: If a light source and an observer at 



distances L, respectively D behind, respectively in front of the texture are perfectly aligned with the 
texture, the observer sees an Einstein ring with a time dependent opening angle 
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This ring opens up at time t = D — eirLD/^L + D), reaches a maximum opening angle /3max = 
\/27reL/(L + D) at t = D and shrinks back to a point at t = D + eTrLD/{L + D). It exists over a 
time span At = l-Ke^^. For reaUstic values of e ~ 10~^ the maximum opening angle can become 
about 1'. 

If source and observer are not perfectly aligned, but deviate by an angle 7 < (/?max from alignment 
the ring reduces to two points which first move apart and then together again within the time At 
and with a maximum separation angle /^max- Unfortunately, even if textures exist, since late time 
textures are so rare, the probability of observing this effect is rather small. There is typically one 
texture which has collapsed after 2; = 4 and for which the photons that have passed it close to the 
time of collapse are visible for us now (see Durrer and Spergel [1991]). The probability that there 
is a quasar behind it with angular separation less than /3max ~ 10~^ is even smaller than for global 
monopoles (see Section 3.4.1). 
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Chapter 5 

Textures in Expanding Space 



We now want to discuss the effect of collapsing textures in a spatially flat, expanding Friedmann 
Universe. To be specific, we assume a universe which is dominated by cold dark matter (CDM) 
which contributes about 95% of the total matter density of the universe. The remaining 5% are 
baryonic. In the next section we write down the system of equations which describes a spherically 
symmetric collapsing texture in this background and the perturbations it induces in the dark matter, 
baryons and photons. We have solved this equations numerically. We then indicate how the single 
texture events can be distributed in space and time to obtain a full sky map of the cosmic microwave 
background. We report briefly on the results of this numerically obtained microwave sky and compare 
it with the COBE data. (A more extensive presentation of this numerical work, which was done in 



collaboration with A. Howard and Z.-H. Zhou will be presented elsewhere [Durrer et al., 1993].) 

Numerical simulations and analytical estimates have already shown that the texture scenario 
leads to early formation of small objects which are likely to reionize the universe as early as z ~ 100 
(see Sect. 4 and Durrer [1990], Spergel et al. [1991] and Gooding et al. [1991]). Photons and baryons 
then are coupled again via Compton scattering of electrons (see Section 3.2.3). They remain so 
until the electrons are too diluted to scatter effectively. This decoupling time is determined by the 
distance into the past, at which the optical depth becomes unity: 

T{tdec) = / Aaarrie = 1 . 

In a flat universe, = 1 one obtains 1 + Zdec ~ 100 ^ niMo ) ^ ' ^'^^^^ denotes the density 
parameter of baryons, and h^o is Hubble's constant in units of 50km/s/Mpc. This last scattering 
surface is of course not as instantaneous as the recombination shell. The decoupling due to dilution 
is a rather gradual process, and the thickness of the last scattering surface is approximately equal to 
the horizon scale at Zdec ~ t dec- 
Due to reionization, perturbations smaller than tdec but larger than the mean free path of the 
photons are exponentially damped by photon diffusion. A rough estimate of this effect is given in 
Section 3.2.3. Results of a numerical calculation of this damping are shown in Fig. 7. 
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5.1 Formalism 



In this section, we present the equations for calculating the response of matter and radiation to the 
collapse of a single spherically symmetric texture. The background metric is given by 

The spherically symmetric ansatz for an texture unwinding at a given time t = tc 'is, like in 
the preceding chapter 

= r/^(sinx sin cos sin X sin sine/?, sin X cos 0, cos x) , (5-1) 
where 9 , (p are the usual polar angles, and xl?'; t) has the properties 

X(r = 0,t<te) =0 
X(r = 0,t>tc) =vr 
x(r = oo,t) = TT . 



In expanding space, the evolution equation (4.5) is replaced by 

dfx + 2{d/a)dtx - dlx - -drX = • (5-2) 

We parametrize the energy momentum tensor of the spherically symmetric texture field, 

T^y = d^4> ■ dycj) - ^g^„ydx(t) ■ d^cj) , (5.3) 
in terms of scalar seed perturbations 



where e = KSirGrf' . Inserting ansatz (5J.) in (^.j), one finds 

fp/f = l/mtX? + {drX? + , (5.4) 

U/f = l/mtxf - \{drX? - ^fj^] , (5.5) 

/-A = 7 ri9tx){drX)dr, (5.6) 



4 



A/, = ka,.xy' - ^1 + 1 rmrxf- - '^f-^ . (5.7) 

4 4 Joo r 

The variables /p//^ and fp/t^ denote the energy density and isotropic pressure of the texture field, 
is the potential of the velocity field and /^r is the potential for anisotropic stresses. Spherically 
symmetric perturbations are of course always of scalar type. Due to the adoption of spherical 
symmetry, we loose, e.g., all information about gravitational waves produced during the collapse. 
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In addition to the texture, we want to describe dark matter, baryons and radiation. Since dark 
matter has zero pressure, we just need a variable describing its density perturbation. 

By Pd we denote the density of dark matter and D is its gauge invariant density perturbation, 
as defined in Chapter 2. The evolution of the dark matter fluctuations is governed by ( p.l35| ) for 
= w = 0: 

D + (d/a)I) - AT.Gpdo'D = e{fp + S/p)//^ , (5.8) 

To describe the baryon- photon system we need two additional variables: The potential for the 
baryon velocity V and the perturbation of the energy integrated photon distribution function M . At 
times t < tdec, the collisionless Boltzmann equation for the photons and the equation of motion for the 
baryons must be modified to take into account scattering. The dominant effect is non~relativistic 
Compton scattering by free electrons. Let us denote the collision integral which is calculated in 
Section 3.2.1 by C{M, V). As in Chapter 3, e is the direction of the photon momentum. One finds 

C = aoTUelDp -M + AeidiV + ^e^.M'^] , (5.9) 

with 

3 



M'^ = —J M{e)eijdQ] ; = e^e,- - (l/3)5i,- 



(We do not have to worry about the position of spatial indices of perturbation variables, they are 
raised and lowered with the Euclidean metric 5ij since fc = in this chapter.) 

The drag force due to Thomson drag of photons on the matter is given by (3.19) 



with Mi = (3/47r) / eiMd^. Including this drag force into the equation of motion for the baryons, 
( ^ ) for w = cl = Q, one obtains dOOp 

IdiV + {a/a)ldiV = di^ - ^^I^{AldiV + Mi) or 

3pb 

lAV + (a/a)lAV = + aaxne—iM - -lAV) . (5.10) 

Pb 3 

For the last equation we have used the zeroth moment of Boltzmann's equation, the continuity 
equation, 

M + {l/3)diM' = . 

In our numerical computations, we have made the simplifying assumption He = for z > Zi and 
ne = for z < Zi for some ionization redshift Zi ~ 200. 

The evolution of the photons is given by the perturbation of Boltzmann's equation, (3.26). 

1 - u2 

M + fidrM + -^^d^M = 

ifidri'^ - ^') + aarUelM -M- 4pdrV + 3(/x2 - 1/3)M2] (5.11) 
1 

where M2{r,t) = - Mip')ip'^ - l/3)dfi' , (5.12) 
2 J -I 
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where fi is the direction cosine of the photon momentum in the direction of r. 

As in Chapter 4, the photon evolution equation is more transparent in characteristic coordinates, 
{t,T,b), where b = r\J\ — [i^ is the impact parameter, and t = t — rfi is the impact time. In these 
variables ( |5.12| ) simplifies: 

dtM{t, T,b) = A^xdri^ + C{t,T,b) , (5-13) 

where C{t,T,b) is the collision integral above, expressed in terms of the new variables. 

In order to write down the perturbed Einstein equations, we have in principle to calculate the 
energy momentum tensor of radiation from A4 . But since we are only interested in late times where 
density perturbations can grow, pd > Pr, we may neglect the contribution of radiation to the density 
perturbation. The potential $ is then determined by the texture and dark matter perturbations 
alone: 

A$ = -e(/p//2 + 3{d/a)fjl) - At^Go" p^D . (5.14) 

On the other hand, since dark matter does not give rise to anisotropic stresses, we have to take 
into account the contribution of radiation to the latter. To calculate the anisotropic stresses of the 
photons we recall the definition of the amplitude of anisotropic stresses, 11: 

The anisotropic contributions to the energy momentum perturbations of the photons are given by 

5Ti-\5Tl5l = ^j{e,e^-\5l)M<m. 
Using these equations and spherical symmetry, one finds 

n" - li'/r = ^ - ^-)Mdp = ^M2(r,t) . (5.15) 

This anisotropy and the anisotropy of the texture field contribute to the sum of the two Bardeen 
potentials 

A(^' + $) = -2A{eU + (4/3)7rGaVrn) . (5.16) 

Let us choose initial conditions that are physically plausible. If the phase transition that produced 
texture occurred in an initially uniform universe, causality requires that there are no geometry 



fluctuations well outside the horizon |Veeraraghavan and Stebbins, 199C]. This implies that initially. 



at ti <C tc, we must require ^' = $ = 0. We want to compensate as much as possible of the initial 
texture fluctuations with an initial dark matter perturbation. Hence, we use as initial conditions for 
the density field, 

D{r, t = U) = - ^^J^2p^ ifp/^' + 3(q/«)/./0 • (5.17) 

This initial condition implies that metric fluctuations are induced by the differences between the 
texture equation of state, and the equation of state of the background matter. This choice yields 
^ = at t = ti, but not ^ = 0. Due to its equation of state the dark matter cannot compensate 
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anisotropic stresses of the texture. We thus must compensate them by an initial photon perturbation. 
For ^' to vanish, we have to require according to ( ^.16 ) 



n = 

This does not lead to a unique initial condition for M, but if we, in addition, require the zeroth and 
first moments of to vanish (which otherwise would interfere with eq. ( ^.17| )) it is reasonable to 
set 

M{r,^^,t = U) = -^^l^y - 1)(A^ - ^f^) . (5.18) 



Together with initial conditions for the texture field xi^^i = U), the requirements ( 5.17| ) and (5.18) 



and the evolution equations ( |5.2D to (|5.16D determine the system which we have solved numerically. 

Figure 9 shows the microwave background fluctuations induced by the collapse of the texture as 
a function of r for small impact parameter in the expanding universe for different times. At t > tc 
it is interesting to see the blueshift at r ~ t of the photons which have fallen into the dark matter 
potential but have not yet climbed out of it again. (During their way out of the dark matter potential 
this blueshift will of course be exactly compensated. This is also visible in the figure.) Figure 10 
shows the microwave background fluctuations induced by texture collapse as a function of impact 
parameter. Note that the temperature fluctuations are induced only for photons that pass within 
the event horizon of the texture. 



5.2 Textures and the Microwave Sky 

In the previous section, we described how the collapse of a single texture produces fluctuations in 
the photon temperature. In this section, we sum the contributions of many textures and describe 
how to construct microwave maps of the night sky. 

Since COBE observations cover the entire celestial sphere, we construct a numerical grid consist- 
ing of 1 square degree patches. These patches are arranged so that they cover equal areas and each 
patch has roughly the same shape. 

Since the texture fluctuations are in the linear regime, we assume that the contributions of each 
texture to the fluctuations at each point in the sky can be added independently. We randomly 
throw down textures everywhere within the event horizon using the texture density distribution 



Bpergel et al., 1991|, 



dn 1 1 , , 

* = 25?- ^'-''^ 

We then follow the collapse of each texture and sum their contributions. We include only textures 
that collapse after recombination. Textures that collapse earlier do not contribute significantly to 
microwave fluctuations on scales accessible to COBE. 

In order to simulate the COBE observations, the map of the night sky is smoothed with a Gaussian 



beam with a FWHM of 7°, the angular resolution of the DMR detector on COBE jSmoot et al., 1991 1. 
After computing and removing the intrinsic dipole contribution (10^^ — 10^'*) and any monopole 
fluctuations, we then compute the microwave quadrupole, the r.m.s. pixel-pixel fluctuations and 
other statistics of the microwave sky. 
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5.3 Results 



Here we present our results for the not reionized universe (zi < z^ec)- We just compare the numerical 
results with the COBE observations on 6 > 10°. In the texture scenario, these large angular scales 
are not affected by reionization. 

The amplitude of the microwave background fluctuations depends upon the scale of symme- 
try breaking associated with the texture. If the symmetry breaking scale is normalized so that 
the scenario can reproduce the amplitude of the galaxy-galaxy correlation function, we have to set 
the value of the dimensionless parameter e = IGirGr]^ to e ~ 5.7 x 10^^6~^, where b, the bias 
factor, is the ratio of the mass- mass correlation function to the galaxy-galaxy correlation function 
Gooding et al., 1991| . Note that the definition of e in this work is 2/7r times the definition used 



in Gooding et al. [1991]. Comparison of the predictions of the texture model with observations of 
clusters [ Bartlett et al., 1993| ] suggest that 6 « 2, a value compatible with hydro dynamical simula- 



tions of texture-seeded galaxy formation | Cen et al., 1991 1. We normalize the microwave background 
fluctuations to this value and present our results in units of eo, where eo = e/2.8 x 10~^. 

We numerically performed 100 realizations of the model. For illustration, a 'COBE map' produced 
from a typical simulation is shown in Fig. 11. Averaging over all realizations, we find an r.m.s. value 
for the quadrupole moment of 

Q = (1.4 ± 1.2) X lO^^eo . 

Since only a handful of textures are the source of most of the large scale fluctuations, the quadrupole 
varies signiflcantly from realization to realization. 

To compare with the COBE result, we have smoothed our calculations over an angular scale of 
10 degrees. The average pixel-to-pixel fluctuations of the smoothed simulations are 

(Ar/r),„,(10°) = (3.8 ± 2.6) X lO-^eo 

The distribution of temperature fluctuations are only mildly non-Gaussian, the skewness of the 
distribution is —4 it 2 and the kurtosis of the distribution is 32 it 29. The errors quoted are statistical 
la deviations of one hundred realizations. One example for the pixel distribution of the fluctuations 
is shown in Fig. 12. 

The results from the COBE differential microwave background radiometers [DMR] place strong 
constraints on CMB fluctuations on scales larger than 10°. In this experiment 

HWright et al., 1992| , gmoot et al., 1992|1 a value of Q « (0.6 ± 0.2) x 10"^ has been found for the 



microwave quadrupole and 

{AT/T)rnU^ = 10°) = (1.1 ± 0.18) X 10"^ 

is the result for the fluctuations at a scale of 10 degrees. The overall spectrum is compatible with a 
Harrison-Zel'dovich spectrum: 

{AT/T){e) « io~5(0/io°)~(°-i^°-^)/2 _ 

Our simulations with bias factor b = 2 are compatible with these results within one sigma. The 
(AT/T) of the simulations is somewhat large but the scatter is considerable. 

We feel that the adoption of spherical symmetry may lead to underestimates of (AT/T) since 
contributions due to gravitational waves and random fluctuations in the scalar field which do not 
give rise to texture (i.e., topological winding number) have been neglected in this approach. If, e.g., 
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{AT/T){6 = 10°) is enhanced by a factor of 2 by these contributions, the COBE measurement is 
about 2.5a below the average texture result. 

The conclusions from these simulations can thus be put as follows: The CMB anisotropies from 
spherically symmetric texture collapse are slightly high but agree within la with the COBE mea- 
surements. This is encouraging. But in the results given above only statistical fluctuations, (i.e., 
cosmic variance) have been taken into account. Due to the uncertainties in modeling the typ- 
ical texture and the approximations inherent in modeling the texture as spherically symmetric, 
these estimates of the microwave background fluctuations are uncertain by at least a factor of ~ 2 
(systematic error), probably underestimating the true induced fluctuations. A full 3d simulation, 
which takes into account also gravitational waves and non-top ological fluctuations of the scalar 
field, is necessary to finally decide on the scenario. Such simulations have now been performed 
[[Bennett and Rhie, 1992| , pen et al., 1993|] . They obt ain results which are higher than those obtained 



in the spherically symmetric approach by a factor 1.5 to 2. 

Smaller scales, « 1° lead to somewhat larger fluctuations and much smaller standard deviations 
(since many textures contribute to them): 

(Ar/r),„,(2°) = (4 ±0.8) X lO-^eo 



If the new measurements which require (Ar/r)(l°) < 1.4 X 10-5 g aier et al., 1992| ] are confirmed. 



Zdec = 100 



we need reionization to damp small scale fluctuations in the texture scenario. 

A rough estimate of the effects of reionization can be obtained by just smoothing each texture 
with a smoothing scale of about the horizon size at Zdec^ 

0.05 

This corresponds to an angular scale of 

Osmooth = tizdec)/to = (1 + z)-^/^ « 5.7° . (5.20) 

If the formation of objects leads to reionization prior to Zdecj this would suppress microwave back- 
ground fluctuations on scales 9 < 9 dec ~ 6°, but would not effect the fluctuations on larger scales 
which are discussed above. 



5.4 Conclusions 

The texture scenario of large scale structure formation has many attractive features. Its galaxy 
galaxy correlation function and the large scale velocity fields agree better with observations than 
in the standard cold dark matter model [ jGooding et al., 1992 1. Also a couple of other statistical 



parameters (Mach number, skewness, kurtosis) are in satisfactory agreement with observations. New 



simulations which used the COBE quadrupole to normalize the fluctuations |Pen et al., 1992] hint 



that, in contrary to earlier results, the scenario may have severe difficulties to reproduce the very 



large scale galaxy clustering observed in the infrared [ Fisher et al., 1993 1, just like the standard CDM 
model. Due to the steep dark matter potential produced with texture, small scale structure can form 
very early. This may lead to early reionization of the universe. Calculations of the microwave 
background anisotropies show that reionization is necessary to reconcile with CMB anisotropies on 
small scales (up to about 5°). Due to the unique signature and relatively large amplitude, the CMB 
anisotropies produced by textures are one of the most hopeful criterion for confirming or ruling out 
this scenario. 
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Another interesting observational test, lensing of background quasars by a foreground texture, is 
very improbable. 

This scenario certainly deserves further work. Especially the investigation of the question if 
reionization as early as z ~ 200 — 100 is possible and a careful analysis of CMB anisotropics for the 
reionized model on angular scales around 1° are important tasks. 
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Appendix A 

The 3 + 1 Formalism 



Since we have made use of 3 + 1 split of general relativity in deriving the linear cosmological per- 
turbation equations in this review, we want to derive the necessary tools in this appendix. A math- 
ematically rigorous overview is given in Choquet-Bruhat and York [1980] and Fischer and Marsden 
[1979] . But explicit calculations are missing there, we shall thus be rather detailed. In this appendix 
we mainly follow the exposition in Durrer and Straumann [1988]. 



A.l Generalities 

We assume that spacetime {M,g) admits a slicing by slices T^, i.e., there is a diffeomorphism 
(j) : A4 — > SxJ, / cR, such that the manifolds = (/>~^(S x {t}) are spacelike and the curves 
(/>~^({x} X I) are timelike. These curves are what we call preferred timelike curves. They define a 
vector field dt, which can be decomposed into normal and parallel components relative to the slicing 
(Figure 13): 

dt = an + l3 . (Al) 

Here n is a unit normal field and /3 is tangent to the slices Ef. The function a is the lapse function 
and f3 is the shift vector field. 

A coordinate system {x*} on S induces natural coordinates on Ai: cl)^^{m,t) has coordinates 
(t, X*) if m £ T, has coordinates x*. The preferred timelike curves have constant spatial coordinates. 
Let us set /3 = /3*0j {di = -^)- From g{n, di) = and (Al) we find 

g{dt, dt) = -a" + (3'(3i , g{dt, di) = f3i . 

In "comoving coordinates" thus 

g = -[a'^ - /3'i3i)dt^ + 2pidx'dt + gijdx'dx^ (A2) 

or 

g = -a'^dt^ + gij{dx' + I3'dt){dx^ + l3Ut) . (A3) 

This shows that the forms dt and dx* -|- are orthogonal. 

The tangent and cotangent spaces of M have two natural decompositions. One is defined by the 
slicing 

Tp{M)=Hp®Vp, (A4) 
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where the "horizontal" space Hp consists of the vectors tangent to the slice through p and the 
"vertical" sub-space is the 1-dimensional space spanned by {dt)p (preferred direction). The dual 



decomposition of (A4) is 



T;{M) = H;ev; , (as) 

with H* = {uj e T*{M) : {oj,dt) = 0} and V* = {uj e T*{M) : {uJ,Hp) = 0}, which is spanned by 
{dt)p. 

The metric defines - through the normal field n - yet another decomposition 
Tp{M) = Hp®H^, (A6) 
where Hp is spanned by n, and dually 

T;{M) = (Vp*)^ © Vp* , (A7) 



Equation (Al) refiects the fact that in general the two directions Vp and H^ do not agree. Dually this 
implies tha t H* and (V*)-'- do not coincide. We have for oj-^ G (^*)"'~ the following decomposition 
relative to ([AS]) 

uj-^ = hor(w-^) + {uj^,l3)dt . (A8) 



The decompositions ( [A 41 ) to ( |A7| ) induce two types of decompositions of arbitrary tensor fields on Ai. 
We call a tensor field horizontal if it vanishes, whenever at least one argument is dt or dt. Relative 
to a comoving coordinate system such a tensor has the form 

S = S'jizf^di^ © • • • (g) dx^^ © • • • © dx^'' . 

This shows that horizontal tensor fields can naturally be identified with families of tensor fields on 
St, or with time-dependent tensor fields on T, ("absolute" space). We denote them with boldface 
letters (except di and dx*). 

As an often occurring example of a decomposition, we consider a horizontal p-form u and its 
exterior derivative du. We have 

du = dcj + dt A dtu> , 

where do? is again horizontal. In comoving coordinates d involves only the dx^ (d = dx^ A di) and 
dfCtJ is the partial time derivative, du and dfCtJ are horizontal and can be interpreted as t-dependent 
forms on E. In this interpretation duj is just the exterior derivative of u. Similarly, other diff'erential 
operators (covariant derivative, Lie derivative, etc) can be decomposed. We use two types of bases of 
vector fields and 1-forms which are adapted to ( |A^ ) and ([A^), respectively ( |A6| ) and ([A7|). Obviously, 
the dual pair {d^} and {dx^} for comoving coordinates {x^} are adapted to (A4) and (A5). On the 
other hand, equations (^^) and (^^) show that the dual pair 

{di,n} and {dx' + p'dt,adt} (A9) 



is adapted to ( [A6D and ([AT]). 

Instead of {di} we also use an orthonormal horizontal basis {ej} {g{ei,ej) = 6ij), together with 
the dual basis instead of {dx*}. Then we have the following two dual pairs, which are constantly 
used: 

{ei,dt} , {-d^dt} (adapted to slicing), (AlO) 
{ei,eo = n} {9^} (adapted to (A6) and (^)) , (All) 
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where the orthonormal tetrad {0^} is given by 

0° = adt , e' = 'd' + p'dt , (A12) 
with here defined hy /3 = . We note also the relation 



eo = n = -{dt - (3'ei 
a 



A. 2 The connection and curvature forms 

We now calculate the connection and curvature forms in the orthonormal basis introduced above. 

A. 2.1 The connection forms 

Prom the first structure equation , 

del" + u%^e'' = () , 

and the definition of the second fundamental form: 

Kij = -m-j , (A13) 

where n denotes the normal field of the slicing, one finds immediately the Gauss' formulas : (remem- 
ber n = eo) 

u\{ej) = u\iej) (A14) 
^°(e,) = -Kij . (A15) 

We define 

dt-d^ = dj-di . (A16) 
Then we can calculate the following quantities : 

a;°(eo) = a-^a\i , (A17) 

^.(eo) = -a-'u;){/3) + - /3j - c) + 4) , (A18) 



and thus, 



K = ^[Lpg-dt9], (A20) 
where the vertical bar | denotes covariant derivation with respect to g . Using the general relation 
dt{detg) = tT{dtg)detg , 
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we find from ( |A20D 

dtvol{g) = (div/3 - atrK)vol(sr) . (A21) 

Let's derive ( [ATtD , ( |AT^ ) and ( [AT9D briefly . 

(i6l° = d{adt) =daAdt = a|ji?* A = a~'^a\i9' A 61° . 

This together with the first structure equation results in ( Af 7| ). (A18) and ( A19| ) are obtained as 
follows: From the first structure equation and ( |A15 ) we conclude 



^e^ ^eo 



-{Kii+ uj\{eo)) . 



We calculate the left hand side of this equation: 

i^. {d& + dtAdte^ + dp^ A dt) 
a-h,^{if^{uj\A'd^) + dt^'-d(3') 



a 



-ir, A 



The symmetric and antisymmetric contribution of the last identity yield the formulas (|A19| ) and 
( [Af8| ) for Kij and uj^j{eo), respectively. 

A. 2. 2 The curvature forms 

We now want to calculate the 3 + 1 split of Rij, Rqj and Gqo- From the second structure equation, 
n'^=dujt, + u;^,Aoji, 



and equations ( |A14| ) to ( |A19D one finds immediately 

n){ek,ei) = ft){ek,ei) + KlKji-KiKjk (Gauss) 

^]{ek,ei) = Kjk\i-Kji\k (Mainardi) . 

We need also the normal components of r?*^- . By the second structure equation we know 

VL'q = -diKijO^) + d{a~^a\ie°) + uj\ A [K^O^ + a-^a\ie'^) . 



(A22) 
(A23) 



A straight forward calculation leads to 



a-^a\ij{e^ A e^) - dK) A 9^ + K){uj\ A9^ - K{e^ A 6^) + uj'jK\ A 6 



which yields ( |A23D and the normal components of 17 g : 

17*o(e,, eo) = a-^a\ + dK){eo) - K>^^.(eo) + KJ - J^K) . 
Prom equations ( A22| ) to ( A24) we can calculate the Ricci tensor with the result: 

Rl3a = ^"Isiea, Ccr) , 



(A24) 
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iioo = ^Aa + a'\dttr{K) - Lptr{Kij)) + trK^ . (A25) 

With help of (^2|) one finds 

Roi = {trK)\i - . (A26) 
For the spatial components we obtain 

Rij = n^i{eo,ej) + n'l{ek,ej) . 
Using ( [A24 ) and ( [A22| ) for the curvature forms leads to 



Rij = Rij + tv{K)K,j - 2Kfj - a-^a\ij - a'^dtKij - LpKij) + KisU;'j{eo) + Kj-.w^eo) .(A27) 
With help of ( |AT^ ), ( |AT9| ) and ( [A20D one can bring ( [A27D into the form 

hoT{Ricci{g)) = Ricci{g) + tr(K)K - 2^^ - a'^{dtK - LpK) - a'^Hess{a) . (A28) 

Using (^25|) and ( [A27| ) we find 

Goo = l/2(-RoO + Ei -Rjj) ('A9Q^ 
= l/2[i? + {iT{K)f - tr(i^2)] . ^"^^^^ 

A. 3 The 3 + 1 split of hydrodynamics 

Calculations similar to those in the last section lead quite rapidly to a 3 + 1 split of hydrodynamics. 
Let us decompose the energy-momentum tensor into horizontal and vertical components: 

T = eeo eo + eo 5 + 5 (g) eo + T . (A30) 

For an ideal fluid with 

T = {p + p)u®u + pg* (A31) 
we find, setting as in special relativity u = 7(eo + v), 7 = (1 — t?^)~^/^, 

e = -i\p + pv'') , (A32) 

S = {p + p)^^v, (A33) 
T = {p + p)-i'^v (^v +pg* . (A34) 
Now we compute V • T for an arbitrary T. From 

Veo(eeo ® eo) = -^^eo(e)eo ® eo + eu;o(eo)ei ® eo + eeo a;o(eo)ei 

and 

Vefe(eeo eo) = Lefc(e)eo (8) eo + eu;o(efc)ej (g> cq + eeo ® t*^o(efc)ei 
we obtain 

V • (eeo ® eo) = Leo(e)eo + eu;o(eo)ej + eu;o(ei)eo ■ 
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In the same manner one finds the other contributions with the result: 

(V • Tf = Leo(e) + ea;*o(ei) + ^° (eo)^^' + S% + u;]ieo)S^ + u;]{ei)T'^ 

Inserting the expressions for the connection forms given in Section Al, leads to the following form 
of the energy equation: 



a 



Similarely one finds 



- V • 5 - 2V(ln a) • 5 + eti{K) + tr(ii: • T) 



and from this we obtain the momentum conservation 

Lp)S = -V(lna)e + 2K -3 + tr(K)5 - a~^V-(aT) . 
This equation is used to derive the vector perturbation equation ( |2.63| ) in Chapter 2. 



a 



(A35) 



(A36) 



A. 4 The 3 + 1 split of Einsteins field equations 

Here we discuss the often used 3+1 split of the gravitational field equations. The calculation of the 
curvature forms relative to the basis ( A12| ) is presented in Section A2. The reader will note that 
Cartan's calculus leads rather quickly to the required results. 

We use the notation introduced in the previous section ( pO|) for the various projections of the 
energy-momentum tensor T into normal and horizontal components. From equations ( [A26| ) , ( |A28 ) 
and ( [A29 ) for the Einstein and Ricci tensors, Einsteins field equations can be written in the form 
(recall that boldface letters always refer to the slices S^): 



R+{tvKf - tvK^ = levrGe , 
V K- V-tT{K) = SttGS , 



(A37) 
(A38) 



5, 



K = Lf^K - Hess(a) + a[Ric{g) - 2K ■ K + (trK)K - 87rG(T - -g{e - trT))] . (A39) 



In addition to ( A37] ), (A38) and ( |A39| ) we have the following relation (Section A2, equation ( |A20| )) 
between g and the second fundamental form K: 



dtg 



(A40) 



-2aK + Li^g . 

Note that this decomposition into constraint equations ( [A37 ), ( [A38 ) and dynamical equations ( [A39D , 
( [A40| ) involves only horizontal quantities and those provides the 3 + 1 split of the gravitational field 
equations. 

In Chapter 2 we also use the following consequence of ( |A39D and (A37) 

dttv{K) = - Aa + L^tT{K) + a[ti{K^) + l/2(e + trT)] . (A41) 

Note that dt and tr do not commute. With (|A40| ) one shows easily 

tiidtK - LpK) = dtti{K) - Lpii{K) + 2ati{K'^) . (A42) 

To derive the perturbation equations for vector perturbations we mainly use ( |A38 ) and ( A39D . 
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A. 5 The 3+1 split of the Liouville operator for a geodesic spray 

In this section we derive a useful form of the Liouvihe operator for a geodesic spray for an arbitrary 
3 + 1 spht. 

We start with some generahties. The metric g of the spacetime manifold M defines a natural 
diffeomorphism between the tangent bundle TAi and the cotangent bundle T*M., which can be 
used to pull back the natural symplectic form on T*Ai. In terms of natural bundle coordinates the 
diffeomorphism is given by ^ {x^,p^ = g^uP^) and those the induced symplectic 2-form on 

TM is 

uj = dx^' A dig^^p'') . (A43) 

The Lagrangian L = i^g^uP^p'^ on Tjv[ defines a Hamiltonian vector field Xg on TA^, determined 
by 

iv OJ = dL . 

^9 



In terms of natural bundle coordinates the geodesic spray Xg is given by 

d_ 

dp^ ' 



Xg=jl^d^-T%p'^p^— (A44) 



where are the Christoffel symbols for {M.^g). (For further details see Stewart, 1971 ].) 

The one-particle phase space for particles of mass m, i.e., the sub-bundle {p G TM : g{p,p) = 
— m^}, is invariant under the geodesic flow and we denote the restriction of Xg to the one-particle 
phase space also by Xg. 

Let / be a distribution function on the one-particle phase space. The Vlasov and Boltzmann 
equations for / involve the Lie derivative Lj^ f. If we consider the spatial components p^, relative 
to an orthonormal tetrad {e^} as independent variables of /, then the Liouville operator Lj^ can 
be written as 

= (A45) 

where w'^ are the connection forms relative to the dual basis {^^}. 

We derive now a more explicit expression of ( A45 ) for an arbitrary 3 + 1 slicing. In order to 
do this, we need the connection forms relative to the basis {0^^} calculated in Section A2 . They 
can be expressed in terms of a, (3,u^j, dj. Using equations ( |A17] ), ( [A15| ), ( |A18 ) and ( A14 ) we find 
(p = p^Ei, _p° = Vp^~+~"?): 

= - Kypi^ + u^^{p)p^^+u;)ie'^)p^p^^ 

= ipOfa-'a\^^ + u^p - a-^(3p°)p^^ - £^(/3/ - c/y |, . 

Here K^j are the components of the second fundamental form of St, for which we also use equation 
([A19|) of Section A2. 

This leads to the following useful 3 + 1 split of the Liouville operator: 

LXgf = i^dt + " [^'^-^^ " i^^^ + (P°)'(lna)l^ - 1^ , (A46) 
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where we have introduced the horizontal tensor field 

H) = a'HP) - c'j) . (A47) 
Equation ( A46| ) is used in Section 2.3. 



A. 6 Glossary 

In this appendix we provide a glossary of the variables used in the text. For most terms we give a 
short explanation and refer to the equation or section where this variable is first used. Usually it is 
defined there. If not, this should be a very common variable found, e.g., in most basic text books on 
general relativity (like the Christoffel symbols, the Riemann tensor and so on). 

A Perturbation of the 00 component of the metric, respectively the lapse function (2. 4), (2. 11), 
Appendix A. 

B Scalar perturbation of the Oi component of the metric, respectively the shift vector (2. 5), (2. 11). 
Bi Vector perturbation of the Oi component of the metric, respectively the shift vector (2. 7), (2. 12). 
Bij Magnetic part of the Weyl tensor (2.27). 

Weyl tensor (2.26,27). 

D*^"^ Gauge invariant density perturbation variable for the matter component a (2.38). 

-Dg"^ Gauge invariant density perturbation variable for the matter component a (2.37). 

-Oi"^ Gauge invariant density perturbation variable for the matter component a (2.36). 

Eij Electrical part of the Weyl tensor (2.26). 

F Gauge dependent perturbation variable for the distribution function, paragraph (2.3.1). 

J^^^^ Gauge invariant perturbation variable for scalar perturbations of the distribution function 
(2.69). 

.F^"^) Gauge invariant perturbation variable for tensor perturbations of the distribution function, 
paragraph (2.3.2). 

!F^^^ Gauge invariant perturbation variable for vector perturbations of the distribution function, 
paragraph (2.3.2). 

G Newtons constant, G = 6.6720 x IQ^^ cm? g^^ sec'"^ . 

G/Mu = RfMu — (l/2)5^i/-R, Einstein tensor. 

H]^ Trace perturbation of the spatial part of the metric. 

Ht Anisotropic scalar perturbation of the spatial part of the metric (2.6,11). 

Hi Anisotropic vector perturbation of the spatial part of the metric (2.8,12). 
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Hij Anisotropic tensor perturbation of the spatial part of the metric (2.9,13). 

Kij Second fundamental form, Section 2.1, Appendix A. 

U Spatial components of the vector field X parametrizing a gauge transformation, Section 2.1, 
2.3. 

Lx Lie derivative w.r.t the vector field X, Section 2.3. 

M Mass used to parametrize the energy momentum tensor of seed perturbations. Section 2.5. 
M. Gauge invariant perturbation variable for the energy integrated photon distribution (2.83). 
M Spacetime manifold. Appendix A, Section 2.3 
Mi = ^J dnciM The first moment of M (3.19). 
Mij = J dQeijM The second moment of M (3.17). 
P^'^ = u^u'^ + (5^ The projection operator onto the 3-space orthogonal to u (2.30). 
Pm The mass bundle, Section 2.3 
R The Ricci scalar. 
R = ^Pm/^Pr Parameter used in paragraph 3.2.3. 

Tl Perturbation of the scalar curvature on the slices of constant time (2.14). 
Rfj^n = i?" The Ricci tensor. 

i?^ ^ai/ The Riemann tensor 

T Temperature of the cosmic background radiation. 

T Temporal component of the vector field X parametrizing a gauge transformation, Section 2.1, 
2.3. 

TM. Tangent space to spacetime, Section 2.3. 

TX Tangent vector field associated to the vector field X, Section 2.3.1 
T^uf^ Scalar contribution to the energy momentum tensor of the seeds (2.117,118,119). 
T^uf"^ Tensor contribution to the energy momentum tensor of the seeds (2.122). 
I/lt^'* Vector contribution to the energy momentum tensor of the seeds (2.120,121). 
V Gauge invariant variable for scalar perturbations of the velocity field (2.35). 
Vi Gauge invariant variable for vector perturbations of the velocity field (2.42). 
X Vector field parametrizing a gauge transformation. Section 2.1, 2.3. 
dpi = gfji Partial derivative (vector field) 
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a Cosmic scale factor, Section 1.1. 

h Impact parameter, Section 3.4, 4.2.2. 

c Speed of light, usually set equal to 1 in this text. 



- \lvl P 1 (ca = \IVa.l Pa) Adiabatic sound speed (of matter component a). Section 1.1. 
Tetrad vector field. Section 2.3, Appendix A. 
/ Distribution function of phase space. Section 2.3. 

/tt Gauge invariant scalar potential parametrizing anisotropic stresses of seeds (2.119). 

/p Gauge invariant perturbation variable parametrizing the energy density of seeds (2.117). 

/p Gauge invariant perturbation variable parametrizing the pressure of seeds (2.119). 

Gauge invariant perturbation variable parametrizing the scalar velocity potential of seeds 
(2.118). 

g^xy Metric of spacetime. Chapter 2. 

h Used to parametrize Bubble's constant Hq = h x 100 gecMpc ' Section 1.1. 
h/j^i, Metric perturbation (2.10). 

h Planck's constant, h = 1.0546 x 10~^^crn^(/,sec^^, usually set equal to 1 in this text. 

k Spatial curvature of a Priedmann universe (1.1). 

k Comoving wave number. Section 3.2, paragraph 3.4.2. 

ks Boltzmann's constant, ks = 1.3807 x 10~^^erg/K, usually set equal to 1 in this text. 

I Length introduced to keep perturbation variables dimensionless, in applications it may be set 
equal to a typical scale of perturbations. Section 2.1. 

Ih = t Comoving size of the horizon, Section 1.3. 
q Redshift corrected energy, paragraph 2.3.1. 
t Conformal time. Section 1.1. 
tx Conformal Thomson mean free path. Section 3.2. 
u Energy velocity field. Section 2.1 

V Scalar velocity potential Section 2.1. 

V Redshift corrected momentum, paragraph 2.3.1. 
Vector peculiar velocity field. Section 2.1. 
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w Enthalpy, Section 1.1. 
w^^^ Gauge invariant vector potential parametrizing anisotropic stresses of seeds (2.121). 

v^f' Gauge invariant vector contribution to the energy flow of seeds (2.120). 
z Cosmological redshift. 

r Gauge invariant entropy perturbation variable, Section 2.1. 
r^;^ Christoffel symbols. Section 2.3. 
A Laplacian. 

A Cosmological constant (1.1). 

n Gauge invariant scalar potential for anisotropic stresses. Section 2.1. 
Ilj Gauge invariant vector potential for anisotropic stresses, Section 2.1. 
Iljj Gauge invariant tensor contribution to anisotropic stresses. Section 2.1. 
S Three dimensional spatial hypersurface, Appendix A. 
$ Gauge invariant scalar potential for geometry perturbations (2.24). 
* Gauge invariant scalar potential for geometry perturbations (2.25). 
ri' Gauge invariant perturbation variable for the fluid vorticity (2.43). 
O''^ Curvature 2-form, Appendix A. 
a Lapse function (2.4), Appendix A. 
(3 Shift vector (2.5), Appendix A. 
7ij Metric of a three space of constant curvature. Section 1.1. 
5 Gauge dependent density perturbation (2.28). 
e* Spatial unit vector (e.g. denoting photon directions). Section 2.3. 
AttGM^ Smallness parameter for the amplitude of seed perturbations, paragraph 2.5.2. 

eijk Three dimensional totally antisymmetric tensor (2.27). 
iCj — 'jij , Section 3.2. 

r] Symmetry breaking scale (4.1). 
01^ Orthonormal tetrad of 1-forms, Appendix A. 

i?' Orthonormal triad of 1-forms on the hypersurfaces of constant time. Appendix A. 
L Isomorphism between the perturbed and unperturbed mass bundles. Section 2.3. 
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i Gauge dependent perturbation variable for the energy integrated photon distribution para- 
graph (2.3.4). 

A Parameter in the scalar field potential (4.1). 

A = (1 + /cr^/4)~^ Conformal factor for the metric of a 3 space of constant curvature, Section 2.3. 

/X Cosine between the photon direction and the radial direction, paragraph 3.2.1. 

tt'^ Orthonormal momentum components Section 2.3. 

TT*- Anisotropic stresses (2.31). 

TTL Gauge dependent pressure perturbation variable (2.31). 
P{a)i P{a) Background energy density of component a. 

a Scalar potential for the shear of the equal time hypersurfaces, extrinsic curvature (2.15). 

(7* Vector potential for the shear of the equal time hypersurfaces, extrinsic curvature (2.17). 
or Thomson cross section, ar = 6.6524 x 10~^^ cm^. 

r Physical time, Section 1.1. 

r Impact time (4.30) 

r Optical depth, Section 5.1. 
T^i, Stress tensor (2.30). 

(f) Scalar field. Chapter 4. 

X Variable parametrizing spherically symmetric scalar field configurations (4.4). 

oj Winding number density of the scalar field. Section 4.1. 
u;J^ Connection forms. Appendix 1. 
ujij Vorticity of a velocity field Section 2.1. 
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FIGURE CAPTIONS 



Fig. 1: The spectrum of the cosmic microwave background radiation as measures by COBE 
(Figure from Mather et al. [1990]). 

Fig. 2: Limits on the CMB anisotropy on different angular scales. The COBE result is the only 
positive detection. All other marks represent 95% confidence upper limits. 

Fig. 3: The linear perturbation spectrum of hot dark matter, for one (2) or three (1) types of 
massive neutrinos. The fluctuations are heavily damped on scales smaller than Xps mpi/m1 
(Figure from Durrer [1989]). 

Fig. 4: Simulations of structure formation with HDM (top right and bottom pictures) compared 
with the corresponding picture from the CfA survey (top left). Triangles are high density regions 
identified as galaxies. One sees that the simulations lead to highly over developed large scale structure 
(Figure from White [1986]). 

Fig. 5: CDM simulations (a) and (b) compared with the CfA survey (c). No striking inconsistencies 
are visible at first sight (Figure from Kolb and Turner [1990]). 

Fig. 6: The angular galaxy galaxy correlation function as measured by the IRAS survey (black 
dots) compared with the predictions from CDM models with h = 0.4 (black line) and h = 0.5 (dotted 
line). The open circles and squares are results from an older analysis of the Lick catalogue (Figure 
from Maddox et al. [1990]). 

Fig. 7 The CMB anisotropy (in units of 10^^) from a spherically symmetric texture collapsing at 
z = 30 (left) and z = 200 (right) respectively as a function of angular separation from the center of 
the texture. This figure is calculated for a universe which reionizes at z = 200. It shows how signals 
from small scale textures are substantially damped and broadened by photon diffusion. 

Fig. 8: The hot spot — cold spot signal of a spherically symmetric collapsing texture in units of 
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e ~ 2.8 X 10^^ . The horizontal variable t = t — r cos 9 denotes the 'impact time' of a photon arriving 
at a distance r from the texture at time t traveling with an angle 6 with respect to the radial direction. 

The hot spot-cold spot is shown for photons with fixed impact parameter b = rsinO ~ O.ltc (tc is the 
time of texture collapse). The signal from the expanding universe att = tc, line (1), and t = l.Stc, line 
(2), is compared with the flat space result (dashed curve). The second peak appearing at t = 1.5tc 
is due to the dark matter potential. 



Fig. 9 As Fig. 8 but for different times with time steps At « 0.25tc- One sees an outgoing wake 
of blue shift at r Ri i. This is caused by photons which have fallen into the dark matter potential 
but have not yet climbed out of it again. This blueshift will of course be completely compensated 
by the redshift these photons will acquire during their way out of the dark matter potential. 

Fig. 10 The CMB perturbation in units of e ~ 2.8 x 10~^ as a function of the impact parameter 
b for fixed r ~ 0.5tc = 10. The signal disappears at an impact parameter b ~ IMc {tc = 20 in the 
units chosen). 

Fig. 11 A simulated COBE map as it might look in a scenario with texture + CDM. The color 
scheme goes from — 4x 10^^ (dark blue) to 1 x 10~^ (deep red). Monopole and dipolc contributions are 
subtracted in this map. A description of how the map is produced (in collaboration with D.N. Spergel 
and A. Howard) is given in the text. 

Fig. 12 The statistical distribution of microwave anisotropies in the texture scenario. The number 
of pixels showing a given anisotropy are counted for one realization of the CMB sky. The distribution 
is slightly non-Gaussian with skewness — 1 and curtosis 3. 

Fig. 13 A 3 + 1 slicing of spacetime M- The family of immersions of S into M is denoted by 
it; it{m) = (f)~'^{m,t). 
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